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A BSTRACT 


The objective of the integrated Digital Flight Control System is to prov.de 
rotational and translational control of the Space Shuttle orbiter in all phases of 
flight, from launch ascent through orbit to entry and touchdown, and during pow- 
ered horizontal flights. The program provides a versatile control system structure 
while maintaining uniform communications with other programs, sensors, and con- 
trol effectors by using an executive routine /functional subroutine format. The pro- 
gram reads all external variables at a single point, copies them into its dedicated 
storage, and then calls the required subroutines in the proper sequence. As a re- 
sult, the flight control program is largely independent of other programs m the 
GN&C computer complex and is equally insensitive to the characteristics of the 

processor configuration. 

Section 1 describes the integrated structure of the control system and the 
DFCS executive routine which embodies that structure. Section 2 deals with the 
input and output. The succeeding sections show the specific estimation and control 
algorithms used in the various mission phases. 
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FORWARD 


This document presents the Shuttle digital flight control design base developed 
by the MSC Guidance and Control Division and the Draper Laboratory Control and 
Flight Dynamics Division. The flight control and attitude control of the Space 
Shuttle from launch through touchdown and for powered horizontal flight is 
addressed. The subject matter reflects different levels of design maturity and 
completeness of documentation for the several different mission phases. It is 
anticipated that this document will b® reissued several times as design and test 
programs are completed. 

The unified Digital Flight Control System (DFCS) structure and philosophy 
presented herein represents a framework for a reasonable evolution of Space 
Shuttle control systems. The Entry, Transition and Cruise flight- control systems 
presented (in the DFCS framework) have undergone detailed tests in the Space 
Shuttle Functional Simulator (SSFS), are considered to be mature designs, and 
represent a demonstration of the benefits in coding economy, structural clarity 
and amenability to straightforward verification of the modular DFCS approach. 
Additions to the Cruise control system to support powered horizontal flight reflect 
efforts still in process; preliminary documentation is presented to outline design 
intent. The framework for two different jet select schemes is presented, on® 
which employs table look up and one which embodies a linear programming approach. 
On-orbit control is incomplete; however, phase-plane designs are presented for 
powered and coasting orbital flight. A TVC flight design using gimballed engines 
is not well developed at this time. Although the boost control is felt to be a mature, 
well tested design, it has not yet been integrated into the DFCS and it is not 
presented in this document. 

Considerable design effort is required to complete the DFCS design. Control 
gains will require change to reflect the new 150K air-frame design. Boost control 
does not yet explicitly reflect a detailed gimballed solid-rocket-motors control 
design philosophy. Transition design is impaired by lack of good trans-sonic 
aerodynamic data. In all cases, the current rigid-body designs must evolve into 
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real-world flexible body flight control designs These designs must he inte- 
grated into the DFCS format and in addition must be re-written in the HA.L higher- 
order- language. 

The flight control system designs will continue to develop with the airtrame 
design. This evolution will be simplified by the common structure and by the 
design approach which has stressed system performance on the low side of accept- 
able performance requirements and gain- stabilized dynamics (avoidance of soph- 
isticated filter concepts). These precepts are expected to permit lower sampling 
rates, moderate computer requirements and some insensitivity to air frame 
structure modeling data. 

The documented DFCS has advanced beyond a paper design: many portions of 
the system are programmed and are running on the SSFS. Test activity is in 
progress to evaluate control/structures and guidance /control interactions: the 
results will lead to further refinement of the definition of structures, control and 
guidance. 
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SECTION 1 


INTEGRATED STRUCTURE: EXECUTIVE ROUTINE 

by 

Peter S. Weissman 


A. Integrated Structure 

The DFCS program consists of an executive routine, several major (functional) 
subroutines, and a number of utility subroutines. The executive routine performs 
certain input processing and "housekeeping” which are required regardless of the 
mission phase and calls the major subroutines in the appropriate order. There is 
a major subroutine for each of the following functions: "pad load" manipulation, 

data read and manipulation, filter ana parameter initialization, filter update, con- 
trol logic, filter pushdown, and parameter estimation. 

This structure of executive routine and functional subroutines gives modu- 
larity to the program; it allows complete flexibility in designing the state estima- 
tion and control algorithms for each of the mission phases, designating new sensors 
or control etieciors, and indeed, redefining the mission phases themselves. At 
the 3aroe time, however, the unified structure which groups together the coding 
for similar functir rs encourages economy both in computer storage, and in engi- 
neering design anc p- iramming effort because it is so easy for different mission 
phases to use comm coding whenever it is appropriate to do so. Another advan 
tage to this structure is that the calis to the functional subroutines indicate ths 
completion of one set of calculations and the beginning of another; hence, they are 
suitable break points for external interrupts, and they facilitate the testing of new 
coding. For example, if modification to the orbit thrust vector control law is con- 
templated while other segments of this phase's computations are unchanged, the 
new coding can be added as a totally separate subroutine, leaving the old coding 
intact. During testing, it is then possible to switch between old and new control 
laws following completion of the first segment of state filtering. This feature 
also could be of utmc mportar.ee for emergency or backup coding that is required 
to be functionally sepai i from the primary coding. 

The unified DFCS structure has additional, managerial advantages that are 
valuable in an enterprise as complex as the space shuttle program. The use in 
different mission phases of the same assumptions about input /output, timing and 


other interfaces is encouraged if not enforced. Furthermore, commonality in con- 
ventions, nomenclature and formats is encouraged, which simplifies and clarifies 
both testing and documentation. 

Not all operations of a control system need to be performed with the same 
frequency. The current version makes allowance for 3 sampling intervals, the 
longer intervals being multiples of the shortest interval. In the current cruising 
flight phase, for example, the fast rate <10 sps) is used for angular control, the 
medium rate (2 sps) is used for velocity control, and the slow rate <1/2 sps) is 
used f_ir altitude control and parameter updating. The medium and slow computa- 
tions are offset, in order to prevent them from occurring on the same control 
cycle. Allowance is made for different sets of sampling intervals in different 
mission phases. 

B, Executive Routine 

An overview of the logic flow through the DFCS is shown in Fig. 1-1. The 
symbols with vertical bars are the functional subroutines; the vertical lines suggest 
the division of the routines, in part, into segments devoted to each of the mission 
phases. The DFCS Executive Routine performs the other functions shown m the 
figure and calls the major subroutines in their proper sequence. 

The operations shown in Fig. 1-1 will be discussed in sequence. The "pad 
load" denotes the values for all parameters that are determined prior to a flight; 
i.e. . that are stored in the flight computer on the launch "pad". The Executive 
Routine reads these values only once - at the beginning of the mission (or the 
beginning of the simulation). 

The following flagwords (external switches) are read on each pass: 

restart required 

mission phase (entry, transition, cruise, etc) 

control mode (manual modes are currently 
defined only in the cruise phase) 

parameter estimation enabled 
a ided - manual - land ing- approach enabled 

The time for the external executive or computer timer to next call the DFCS 
is then updated. 

The phase indication is automatically changed at this point if the conditions 
warrent. Then the Executive Routine sets a flag, ISTART, calling for a DFCS 
initialization if it is the first pass, FLAG 1 is set, or the phase or mode has 
changed since the last pass. If an initialization pass is called for, the Executive 
Routine records the new mission phase and control mode, establishes the sampling 


FLAG 1 
FLAG 2 
FLAG 3 

FLAG 4 
FLAG 5 
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intervals for the new phase, and calls the major subroutine which transters the set of 
pad-loaded constants required for the current phase from the pad-load locations into 
working storage. 

Regardless of whether or not it is an initialisation pass, the Executive Routine 
calls the Input Interface Routine, which reads all the dynamic input equtred by the 
current phase and mode : sensor measurements, navi ration estimates, guidance 

commands, and manual summands. All the dynamic data is read at the same point 
to ensure that aU the values are registered with respect to the other values. The 
parameters are differenced, rotated into other coordinate systems and otherwise 
manipulated so as to be in the proper form for use in the other routines. During 
the appropriate phases, an Air Data Routine is called which manipulates some 
aerodynamic parameters. These two routines are described in Section 2. 

At this ooint, the logic flow is detoured in an initialization pass to go through 
the Filter and Parameter Initialization Routine and the Parameter Estimation Routine. 
The former assigns values to those variables that need to be initialized when the phase 
(or mode) is entered. The latter determines the current values of parameters that 
Wl H be regularly updated in the course of this phase (or mode). 

The flow then returns to the normal path the Filter Lpdate Routine, Control 
Routine, and Filter Pushdown Routine are called in sequence. These three sub- 
routines comprise the core ot the control system. Hfter they are executed, the 
Parameter Estimation Routine i* called; it is not called at this point in an mitializa 
tion pass because it was execute! earlier in the pass. The state filters are divided 
into two parts m order to minimize the amount of computation between data mput 
and control command. The first part is limited to updating the filter estimate 
according to the mod current measurement. The second part, which occurs after 
the control law, is used to pushdown data bn the case of classical digital filters) or 
to p. opagaie the state un the ca.se ui disc rote-time modern control estimation). 

The Control Routine scales and combines the error signals from the filter routine to 
form control signals for the aerodynamic sunaces and reaction control thrusters. 

The Parameter Estimation Routine updates the gams required for estimation and 
control. The functions performed in this routine may include scheduled or adaptive 
gam adjustment. In the case of the manual cruising flight modes, this section of the 
program is also used to generate a reference trajectory for automatic speed brake 
control and aided-glide-^iope following. Much of the Parameter Estimation Routine 
is executed only at the slow sampling rate. 

The algorithms incorporated in these routines for the Cruise, transition and 
Entrv phases are presented in Sections 3, 4 and 5. A preliminary design for con- 
trol in the On-Orbit phase is presented in Section 6. 
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After these routines have been executed, the Executive Routine increments the 
indices that indicate the "medium-" and "slow- sampling- rate" DFCS passes. " •: 
the end of an initialization pass, the [START flagword is reset and the "slows. ,.p- 
ling-rate" index is given an extra increment. This extra increment will smooth 
the load on the computer by offsetting future "slow" passes from "medium" passes, 
provided that one sampling rate is an integer multiple of the other. The executive 
routine then returns control to the external executive. 

Table 1-1 describes the mission phases for which the DFCS currently makes 


provision. 

Table 1-2 and 1-3 define the symbols and the conventions which are used in 
presenting the DFCS equations in this document. 
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Table 1-1 DFCS Mission Phases 


2 


3 


4 


5 

6 

7 

8 


Rollout 


Transition from aerodynamic control to steering and braking 
with landing gear, after touchdown; not currently implemented. 

Cruise 


Approach and landing following entry and transition phases 
(currently below Mach 0. 9); being expanded to include 
powered flight (take-off through landing). 

T ransition 

Transition from high to low angle of attack; utilizes both 
A CPS and ACS; currently from Mach 5 to Mach 0. 9. 

Entry- 


From after the de-orbit maneuver to transition phase; from 
total reliance on ACPS to blended control with ACS. 

Orbit, 

Orbit, 

T VC 
RCS 

To be combined into a single phase; not currently implemented 
studies in progress. 

Insertion 

Insertion into orbit following the boost phase; noi currently 
implemented; studies in progress. 

Boost 


Flight in the atmosphere with the Solid Rocket Motors; not 
currently implemented; studies in progress. 
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Table 1-2: Nomenclature 


Phase + 

Mnemonic 

Svmbol in Text 

in Program 


2 

ACCN 

a 

4 

A DOT 

at 

4 

AEOLD 

" out old 

2.3 

Aid) 

^in 

o 

Aid) 

a. 

l sum 

3 

ALFINT 

f<* 

3,4 

ALFLM 

a . 
out lim 

2, 3.4 

ALPH 

Q 

2,3.4 

ALPHG 

*G 

2,4 

AMEAN 

a u 

2.3,4 

AO(l) 

°o\xt 

2 

A 0(2) 

a 

o 

sum 

2 

AT RIM 

6a tr 

4 

AX, AY, AZ 

U x ,Uy, t 

4 

Al, A 2, A3.A4 

a a^, 

3.4 

ROOT 

3 

4 

BDTM 

3 M 

2.3.4 

BETA 

8 

2.3,4 

BETAG 

3 G 

2 

BI(1) 

8 in 

2 

BT(2) 

3 i 

sum 

3 

BOd) 

3 out 


* Phases are defined in Table 1-1. 


Definition 

Input normal acceleration (measurement) 
Derivative of angle-of-attack WRT time 
Saved value of angle-of-attack filter output 

Angle- of-attack filter input 
Angle-of-attack filter past sum 

Integral of angle-of-attack filter output 
Limited value of angle-of-attack filter outpu* 

Angle of attack (measurement) 

Angle of attack (guidance) 

Mean angle of attack 
Angle-ot-attack filter output 
Angle-of-attack filter past sum 

'Aileron' trim discrete 

jet acceleration in roll, pitch, yaw 

Gains for ’aileron' control 

Derivative of sideslip WRT time 
Rate limit in roll jet logic 

Sideslip angle (measurement) 

Sideslip angle (guidance) 

Sideslip-angle filter input 

Sideslip-angle filter past sum 

Sideslip-angle filter output 
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Table 1-2: Nomenclature (Cont) 


Phase 

Mnemonic 
in Program 

Svmho] in Text 

2 

BO<2> 

3 o 

sum 

4 

CAT 

cos a 

T 

3 

COSA(l) 

cos °M 



* 

3 

COVARtt) 

cv 

2,3,4 

CRT B 

c b 

r 

2,3, 4 

DAC 

6a c 

4 

DACMAX. 

DACMIN 

. max , min 

6a c # 6a c 

2,3 

DAG 

6a G 

3 

DALFIN 

A / i* 

2 

DAM 

6a 

m 

4 

DAT 


4 

DAI, DA2 

^ ct2 

4 

DB1, DB2 

6 01 a 32 

2,3,4 

DEC 

6e c 

2 

DECH 

6 eh 

c 

4 

DECMAX, 

DECMIN 

r max , min 
6e c , oe^ 

2 

DECU 

6eu c 

2,3 

DEG 

6e G 

2,3,4 

DELC 

6L C 

2 

DEM 


2, 3, 4 

DTRC 

6R C 

4 

DET 

6e t 

2, 3 

DMIN(I) 

6 (min) 

2, 3 

DMAX(I) 

6(max) 


Defimuon 

Sides lip 'angle filter past sain 
Cos (A MEAN) 

Cos (AMEAN) approximation 

lateral covariance matrix (upper triangle) 
Reference -to- body direction cosine matrix 
’Aileron’ angle command output 

Current ‘aileron’ command limits 

'Aileron* angle command (guidance) 

Increment in ALFTNT 

’Aileron' angle command (manual) 

'Aileron’ trim estimate 

Deadbands for pitch jet logic 

Deadbands for roll jet logic 

’Elevator’ angle command output 

‘Elevator’ command due to altitude feedback 

Current 'elevator’ command limits 
’Elevator’ command due to velocity feedback 
'Elevator' angle command (guidance) 

Left eleven angle (command) 

’Elevator' angle command (manual) 

Right elevon angle (command) 

•Elevator* trim estimate 
Aero, surface deflection limits 
Aero, surface deflection limits 
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Table 1-2: Nomenclature (Cont) 


Phase 

3 

3 

Mnemonic 
in Program 

DOMEGA(I) 

DPHINT 

Symbol in Text 

A „ 

roll 

A ^ u 

pitch 

A u> 

yaw 

A/* 

4 

DP1STR, 

DP2STR 

* * 

a dpl’ q dp2 

2, 3 

DRC 

6r C 

2, 3 

DRG 

6r 

u 

2 

DRM 

6r 

m 

4 

DR1, DR2 

6 01 ’ 6 <t> 2 

2, 3, 4 

DSBC 

6sb c 

2, 3 

DSBG 

6sb G 

2 

DSBM 

6sb 

m 

2, 3, 4 

DYNP 

q dp 

4 

DYNP,, DYNP2, o dp2 . 

q dp3 

2 

ELBIAS 

6e B 

2 

ETR1M 

6e tr 

- 

FLAG1 

FLAG 1 


FLAG2 

FLAG 2 


Definition 

Rate change commanded by jet law (roll, 
pitch, yaw) 


Increment in PHOINT 

Dynamic pressure for full ’elevator' 

(’aileron*) use 

Rudder angle command output 
Rudder angle command (guidance) 
Rudder angle command (manual) 


Deadbands for yaw jet logic 
Speed brake angle command output 

Speed brake angle command (guidance) 

Speed brake angle command (manual) 

Dynamic pressure estimate 

Dynamic pressure switch levels for ’elevator’ 

•aileron 1 , roll jet mode III 

'Elevator’ bias to trim speed braice bias 

'Elevator* trim discrete 

New start discrete (external) : 

0 - normal cycle 

1 = new start 

Flight phase flag (external) : 

1 = rollout phase 

2 = cruising phase 

3 = transition phase 

4 = entry phase 

5 - orbit TVC phase 

6 = orbit RCS phase 

7 = insertion TVC phase 

8 = booster TVC phase 
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Table 1-2: Nomenclature (Cont) 


Phase 

Mnemonic 
in Program 

Symbol in Text 

Definition 


FLAG 3 


FLAG 3 

Mode flag (external) : 

1 - direct 

2 = rate command 
3 = RCAH 

4 - auto 

- 

FLAG4 


FLAG 4 

Parameter estimation flag (external) * 

0 - skip 

1 = do 

- 

FLAG 5 


FLAG 5 

Manual mode auto, speed brake flag (external): 

0 * no auto, speed brake /Mach trim 

1 = auto. speed brake /Mach trim 

2,3 

G(I) 


G I 

Filter gains for current flight phase /mode 
(in Phase 3, also includes state extrapola- 
tion constants and long, process noise com- 
ponents) 

2.4 

GAIR 


g a 

Filter gain for air data (=TF/TS) 

2,4 

H 


H 

Altitude above runway (measurement) 

2 

HG 


H G 

Altitude above runway (guidance) 

2 

Hid) 


H ln 

Altitude-above -runway filter input 

2 

HI(2) 


H. 

'sum 

Altitude-above-runway filter past sum 

2 

HO(l) 


H out 

Altitude-above-runway (filter output) 

2 

HO(2) 


H o 

sum 

Ahiiude-above-runway filter past sum 


INGUID 


INGUID 

Flag specifies guidance command set 
neg = commands not updated this pass 
+ /- 1 = ROLLG, ALPHG specified 
+ /- 2 = ROLLG, QG specified 


INSENS 


INSENS 

Flag specifies source of measurements 

0 = provided by simulator (hybrid sim.) 

1 = sensors or environment (SSFS) 

2 = G&N and sensor J or environment (SSFS) 

2, 3 

K(I) 


*1 

Control gains for current flight phase /mode 

4 

KA, KE 


K a' K e 

’Aileron* (’elevator*) deflection limit parameter 

4 

4 

KALFT, KPHIT 

LAMDAA, 

LAM DAE 

x a' X e 

Gain for alpha (phi) integrator for trim 

Switch set to 1 if pitch (longitudinal) ACPS 
jets are needed 
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Table 1-2: Nomenclature (Cont) 


Phase 

Mnemonic 
in Program 

Symbol in Text 

4 

MA, ME 

m a* m e 

2,3,4 

MACH 

m 

4 

MU 

u 

2 

NI(l) 

a in 

2 

NIM 

a i 

sum 

2 

NO(l) 

a out 

2 

NO(2) 

a o 

sum 

2,3,4 

P 

P 

2,3 

PG 

PG 

2,3,4 

PHI 

* 

2,3,4 

PHIG 

*c 

2 

PHIP 


2,3 

PHJ(l) 


2 

PHJ(2) 


2, 3 

PHCXl) 

6 

v out 

2 

PH 0(2) 

*o 

sum 

3 

PHO(2) 

^out 

3 

PHOINT 

I* 

3,4 

PHVLIM 

d> 

v lim 

2, 3 

PI(1) 

Pin 

2 

PI(2) 

Pi sum 

3 

PN(1) 

PN 

2,3 

PO(l) 

P out 

2 

PO(2) 

°sum 


Definition 

Slope of 'aileron' ('elevator') deflection 
limit parameter curve 

Mach number 

Intermediate variable in roll jet logic 
Normal acceleration (filter input) 

Normal acceleration filter past sum 

Normal acceleration (filter output) 

Normal acceleration filter past sum 

Roll rate (measurement) 

Roll rate (guidance) 

Roll attitude (measurement) 

Roll attitude (guidance) 

Roll angle in steady turn 

Roll angle (filter input) 

Roll angle filter past sum 

Roll angle (filter output) 

Roll angle filter past sum 

Roll angle filter output extrapolated 

Integral of roll filter output 
Limited value of roll filter output 

Roll rate (filter input) 

Roll ra + e filter past sum 

Noise matrix (upper triangle); added 
to lateral covariance 

Roll rate filter output 
Roll rate filter past sum 
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Table 1-2: Nomenclature (Cent) 


Phase 

Mnemonic 
in Program 

Symbol in Text 

3 

PO(2) 

p out 

2 

PRI 

p 'i„ 

2 

PRO 

Pr out 

2 

PRI 

Pr. 

x sum 

2 

PR 2 

p r 

o 

sum 

2,3 

PSI 

4* 

2,3 

PSIG 

*G 

2,3 

PSJ(!) 

* in 

2 

PSJ(2) 

sum 

2.3 

PSO(l) 

^ out 

2 

PSO(2) 

* o 

sum 

3 

PS 0(2) 

<aut 

2 

FT RIM 

p tr 

2 

PTRl 

p tr 

sum 

2,3,4 

0 

q 

2 

QB 

q B 

2,3 

QG 

q G 

2, 3 

QI(1) 

q in 

2 

01(2! 

^sum 

3 

QUMM 

qout lim 

2, 3 

00(1) 

^out 

2 

QCX2) 

% 

sum 

3 

QO(2) 

q out 

2 

QT 

q t 


Definition 

Roll rate filter output (extrapolated) 

Roll- rate-to- rudder filter input 

Roll- rate- to- rudder filter output 

Roll- rate-to- rudder filter past sum 

Roll-rate-to-rudder filter past sum 

Yaw attitude (measurement) 

Yaw attitude (guidance) 

Yaw attitude (filter input) 

Yaw attitude filter past sum 

Yaw attitude (filter output) 

Yaw attitude filter past sum 

Yaw attitude filter output (extrapolated) 
Roll axis trim command 
Trim filter (temp) 

Pitch rate (measurement) 

Pitch rate comp, in steady turn 

Pitch rate (guidance) 

Pitch rate (filter input) 

Pitch rate filter past sum 

Limited value of pitch rate filter output 

Pitch rate (filter output) 

Pitch rate filter past sum 

Pitch rate filter output (extrapolated) 
Elevator trim bias 
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T*bie 1-2: Nomenclature (Coat) 


Phase 


2.3.4 

4 


Mnemonic 

iri~Program 

QT RIM 

QTR! 

QVAR 

R 

RDOT 


Symbol in Text 


4 tr 


4 tr 


sum 


Definition 

Pitch axis trim command 
Trim filter (temp) 

Pitch rate measurement variance 

Yaw rate (measurement) 
Derivative of aerodynamic roll 


2.3 

RG 

r G 

Yaw rate (guidance) 

2 

RGO 

R GO 

Range to go 

2 

RGOVEC 

-GO 

Input vector range to go 

2. 3 

Rl(l) 

r in 

Yaw rate (filter input) 

2 

RI(2) 

r i 

sum 

Yaw rate filter past sum 

2, 3 

RO(l) 

r out 

Yaw rate (filter output) 

2 

RO(2) 

r 

°sum 

Yaw rate filter past sum 

3 

RCX2) 

r out 

Yaw rate filter output (extrapolated) 

2,3.4 

ROLL 

4> v or <t>* 

Computed roll angle 

2.3,4 

ROLLG 

<t> or <t > . * 

V G b G 

Guidance commanded roll angle 

2,3,4 

ROLLI(l) 

* v or *b* 
in in 

Roll filter input 

2.3,4 

ROLLO(l) 

^ v out 

Roll filter output 

2,3,4 

RR 

R rel 

Position vector in reference coordinates 

2 

RTRIM 

r tr 

Yaw axis trim command 

2 

RTR1 

tr sum 

Trim filter (temp) 

2 

RUTRIM 

6r tr 

Rudder trim discrete 

2. 3 

RWASH 

r wash 

Yaw rate (washout filter output) 

2, 3 

RW1 

r _ 
sum 

Yaw rate washout filter past sum 

4 

R1 

r l 

Angle - of -a ttac k error limits 

Ambiguity 
to the air- 

is resolved by the discrete INGUID. 

When INGUID = 1, the roll is referenced 

-relative velocity vector, and the subscript v applies. When LNULiU - z, me 

roll is referenced to the body 

x-axis, and the subscript b applies. 
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Table 1-2: Nomenclature tConO 


Phase 

Mnemonic 
in Program 

Symbol m Text 

4 

R2, R:i 

* 2* r :? 


R4, R3, R6, 

R7 1 4» r :V r g» 

4 

R8 

r 8 

4 

R'J 

r 9 

4 

RIO, Rll 

r 10* r ii 

3 

SAG 

S' 

out h-l 

4 

SAT 


2 

3BH 

6sbh^, 

2 

SB1AS 


2 

SBU 

6sbu^ 

4 

SIGMA A, 

B,R ° a -%- a 

3 

SINA( l ) 

sina M 

3 

SPH 

6 

V ° ut n-1 

2,3 

T AN PHI 

tan^p 

2,4 

TEMPER 

T A 

2,4 

TEMPRK 

TKq 

- 

TF 

t f 

2,3 

THETA 

9 

2,3 

THETAG 

9 g 

2,3 

THI(l) 

9. 

m 

2 

THI(2) 

9. 

1 sum 

2 

THO(l) 

^ cut 


Definit ion 

y s i c a - elevator' 'ow^r fuot K*r) t ini it 
Parameters for blending elevator . pitch jets 
Roll angle error limits 
Physical 'aileron' deflection limits 
Parameters for blending aileron/ lateral jets 
Angle- of attack filter output on previous pass 

Sin (A MEAN) 

Speed brake command due to altitude feedback 
Speed brake bias 

Speed brake command due to velocity feedback 
Gains for jet logic - pitch, roll, yaw 

Sin (AMEAN) approximation 

Roll angle filter outout on previous pass 

Tan (roll angle in steady turn) 

Air temperature. Deg F and dimensionless 
Input air temperature in Deg Kelvin 
Fast sampling interval 

Pitch angle (measurement) 

Pitch angle (guidance) 

Pitch-attitude filter input 

Pitch-attitude filter past sum 

Pitch-attitude filter output 
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Table- 1-2: Nomenclature (Cent) 


Phase 

Mnemonic 

yrobol in Text 

Definition 

2 

in I'r ogram — 

THO{ 2) 

°surn 

Pitch-attitucie filter* past sum 

3 

THO(2) 

f)' 

out 

Pitch- attitude filter output (extrapolated) 

3 

THVAR 

2 

CT 8 

Pitch angle measurement variance 

- 

TS 

T S 

SI ow sampling interval 

2, 3, 4 

U 

u 

Earth- relative velocity magnitude (measure- 
ment) 

2,3 

UG 

U G 

Earth- relative velocity magnitude (guidance) 

2,3 

UK 1) 

U in 

Earth- relative velocity magnitude filter 
in put 

2 

UI{2) 

Ui ,„ 

Earth- relative velocity magnitude filter past 

sum 

4 

UJ(1) 

u , u , u 
x y z 

Signum of required jet torque for each axis 

4 

UL(I) 

VL 

Vector in longitudinal plane normal to velocity 

£4 

UMEAN 

U M 

Mean velocity 

2,3 

UO(l) 

U out 

Earth- relative velocity magnitude filter- 
output 

2 

U 0( 2) 

L o 

sum 

Earth- relative velocity magnitude filter past 
sum 

2,3,4 

UXV, L YV, L'ZV 

uv ,uv , uv 

— x — y — z 

Unit vectors along the stability axes, in 
reference coordinates 

2, 3 

U1 , U2 

U i- l 2 

Gains for roll jet logic 

3 

VARMES(I) 

<7 a 2 » a 2 

t|f ' r ' 0 p 

Measurement variances of lateral state 
components 

2,3,4 

V BODY 

- body 

Earth- rel velocity in body coordinates 

2,3,4 

VR 

-rel 

Earth- rel velocity in reference coordinates 

4 

WD 

w d 

Gain for elevator control 

4 

W1,W2 


Gains for aileron control logic 

4 

XII - XI 6 


Scheduled control gains 

4 

ZETAD 

f 

’ d 

Gain for elevator control 

4 

ZETA1, 2 

r r 

* 1 • * 2 

Gains for aileron control logic 
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Table 1-3: Symbol Conventions 


S ymbo l 


Definition 


1 ) 

< >' 
v or v 

*1 

c a 


unit ( v } 


tan * ( x ) 


. -1 , , 
sin ( x ) 


< >4 


in 


( >4 


sum 


< >out 

< )o 


sum 


( V 


( ), 


measured quantity 
extrapolated quantity 
vector v 
matrix M 

direction cosine matrix relating frame b to frame a 

magnitude of vector v 
absolute value of scalar a 
vector cross (outer) product 
vector dot (inner) product 
function which unitizes vector v 

arctangent of x 

arcsine of x 
filter input 

filter input past sum 

filter output 

filter output past sum 

guidance command (DFCS input) 

DFCS command (DFCS output) 
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SECTION 2 


DECS INPUT AND OUTPUT 


by 


j. Ed'* in Jones, Rotor 5. WeUsman, Craig C. Work, 
Edward T. Kubiak (MSC/EG2) 


Inputs from G&N and Sensors: Input Interfa ce Routine 


The Input Interface Routine receives DFCS inputs and computes quantities used in 
the state and parameter estimation routines from information supplied by guidance and 
navigation (G&N), information taken from various environment (ENV) modules (in lieu of 
information which would normally be supplied by sensors), and information generated 
within +he DFCS. Inputs to the routine are listed in Table 2-1, aiorq with the associated 
source module. All angular inputs are in radians. All lengths received from G&N are in 
feet, and all lengths received from ENV are in meters. Ail input vectors are converted 
to inertial coordinates except B co . which is in runway coordinates. Outputs of the 
routine are similarly listed in Table 2-2. All angular outputs are in radians and all 
lengths are in feet. Interface -routine operations are governed by the discretes 
INGL'ID, INSENS, ISTART, and MODE in the manner summarized in Table 2-3. 


A logic-flow diagram of the interface routine is included as Fig. 2-1. The logic 
is divided roughly into three sections: data input, command processing, and filter- 
input computation. A discussion of this logic follows. 

The routine is called by the executive routine. If ISTART = 1, signifying an 
initialization pass, guidance-command variables which are not currently used are 
set to zero. 

The data input section of the routine is now entered unless INSENS-0, in which 
case necessary inputs are assumed to be provided directly by the simulator and the 
data input section is bypassed. Otherwise, measurements of the body rates, temper* 
erature in degrees Kelvin, and the reference -to-body direction cosine matrix are 
read from the appropriate environment module. The dimensionless temperature, 

T , is computed from the temperature measurement, TK^ . Next, the air-relative 
velocity, position, altitude, and velocity magnitude are obtained either directly from 
the environment or from G&N, based on INSENS. When INSENS* 2 angle-of- 
attack, a , and angle-of-sideslip, $ , are also read from the environment. When 
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Table 2-1: Inputs to the DFCS Input Interface Routine 

(Source ambiguities are resolved via the 
discrete LNSENS as per Table 2-3 ) 


Svmboi 

DFCS 

Desc ription 

Source 

Mnemonic 




1ST ART 

initialization cycle discrete 

DFCS 


IN SENS 

sensor specification discrete 

External 


INGUID 

guidance- command specification disci et^ 

G&N 

a 

ALPH 

angle- of-attack (measurement) 

G&N or ENV 


BETA 

sideslip angle (measurement) 

G&N or ENV 


ROLL 

roll angle about body x-axis (measurement) 


P 

P 

roll rate (measurement) 

ENV 

q 

Q 

pitch rate (measurement) 

ENV 

7 

R 

yaw rate (measurement) 

ENV 

u 

U 

earth- relative velocity magnitude 
(measurement) 

G&N or ENV 

H 

H 

altitude above runway (measurement) 

G&N or ENV 

tk a 

TEMPRK 

input air temperature in deg Kelvin 

ENV 

-GO 

RGOVEC 

input vector range to go 

G&N 

*rel 

VR 

earth- relative velocity in reference 
coordinates 

G&N or ENV 

-rel 

RR 

position vector in reference coordinates 

G&N or ENV 

C b 

CRTB 

reference-to-body direction cosine 

ENV 

r 


matrix 

(constant , 

a 

ACC'S 

input normal acceleration (measurement) 

not used) 

°vG 1 

ROLLG 

guidance -commanded roll angle 

G&N 

°bG 1 




°G 

ALPHG 

guidance- commanded angle- of -attack 

G&N 

q G 

QG 

guidance- commanded pitch rate (guidance) 

G&N 

LV 

UG 

earth- relative velocity magnitude 

DFCS 

G 


(guidance) 


h g 

HG 

altitude above runway (guidance) 

DFCS 


PHASE 

flight phase flag (internal) 

DFCS 


MODE 

mode flag (within phase) (internal) 

DFCS 

T 
1 F 

TF 

fast sampling interval 

DFCS 

q dp 

DYNP* 

dynamic pressure (measurement) 

ENV 

6 sb G 

DSBG* 

guidance- commanded speed-brake 
setting 

G&N 


''Not manipulated by Input Interface Routine. 
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Table 2-2: Outputs from Input Interface Routine 


Symbol 

DFCS 

Dv* 

c ription 

DVCS Subroutine 
Which Uses the Output 

Mnemonic 





a 

ALPH 

angle-of-attack (measurement) 

Initialization, 
Air Data 

8 

BETA 

sideslip angle (measurement) 

Filter Update, Control 

P 

P 

roll rate (measurement) 

Initialization. Filter 
Update. Control 

q 

Q 

pitch rate (measurement) 

Initialization. Filter 
Update. Control 

r 

R 

yaw rate (measurement) 

Initialization, Filter 
Update, Control 

u 

U 

earth - relative velocity 
magnitude (measurement) 

Initialization. Parameter 
Estimation, Air Data 

H 

H 

altitude above runway 

Parameter Estimation. 

(measurement) 

Air Data 

P n 

PG 

roll rate (guidance) 

Initialization 

q G 

QG 

pitch rate (guidance) 

Initialization 

r G 

RG 

yaw rate (guidance) 

Initialization 

T. 

TEMPER 

air temperature in deg F 

Air Data 

A 

(dimensionless) 


r go 

RGO 

range-to-go 

Initialization, 
Parameter Estimation 

*p 

PH IP 

roll angle in steady turn 

Filter Update 

q isum 

QI(2) 

pitch rate (filter input) 
past sum 

Filter Update 

q 

QO<2) 

pitch rate (filter output) 

Filter Update 

M o smm 

past sum 


or. 

m 

Al(l) 

angle -of- attack filter input 

Initialization, Filter 
Update, Filter Pushdown 


BI(1) 

sideslip angle filter input 

Initialization, Filter 
Update. Filter Pushdown 

^vin | 

ROLLI(l) 

(about Y re i> 

roll filter input 

Initialization, Filter 
Update, Filter Pushdown 

0, . ) 


(about body x-axis) 

y b in 

PI(1) 

roll rate (filter input) 

Initialization, Filter 

*in 

Update, Filter Pushdown 


QIG) 

pitch rate (filter input) 

Initialization, Filter 
Update, Filter Pushdown 
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Table 2-2: Outputs from Input Interface Routine (Cont) 


Symbol 

DFCS 

Mnemonic 

hi m 

Description 

DFCS Subroutine 
Which Uses tne Output 

r 

in 

yaw rate filter input 

Initialization. Filter 
Update. Filter Pushdown 

t * 

IT (1) 

ea rth- relat ive ve kx'ity 

Initialization. Filter 

in 

magnitude filter input 

Update, Filter Pushdown 

H r 

HI (1) 

altitude above runway 

Initialization Filter 

(filter input) 

Update. Filter Pushdown 

a 

NI (1) 

normal acceleration 

(unused) 

in 


(filter input) 



DYNP 

dynamic pressure 
(measurement) 

Air Data 

6sb~ 

DSBG 

guidance - command speed 

Control 

u 


brake setting 
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Table 2-3: Discretes Governing Input - Interface Operations 


M n e monk’ Va luc 


Action 


INGUID 


I 


negative 


Process pitch- rate and roll (about body x-axis) 
guidance commands this pass. 

Process angle- of- attack and roll (about relative 
velocity vector) guidance commands this pass. 
Do not process guidance commands this pass. 


INSENS 2 
1 
0 


When there is a choice* obtain routine inputs 
from the environment. 

When there is a choice, obtain routine inputs 

from GfcN. ^ 

All routine inputs provided by the simulator. 


1ST ART 1 
0 


Initialize internal variables this pass. 

Do not initialize internal variables this pass. 


MODE 4 

<4 


Compute automatic -mode filter inputs. 

Do not compute automatic- mode filter inputs. 


* See Table 2-1, column headed "Source”. 

option is provided for compatibility with the Charles Stark Draper 
Laboratory Hybrid Simulator and similar systems. 


2-5 



Fig. 2-1 Logic Flow Diagram for DFCS Input Interface 
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MODE - 4 
(auto) 



Fig, 2-1 Logic Flow Diagram for DFCS Input Interface (Cont) 



* 




¥ 
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0 



Fig. 2-1 Logic Flow Diagram for DFCS Input Interface (Cont) 
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INSENS* 1. a and JT are computed using the air- relative velocity and the direction 
cosine matrix. C*. Additionally, the range-to-go magnitude (Rqq) is computed 
when INSENS* 1. 

If the DFCS is not in the automatic mode, the remainder of the data input 
section, the command processing section, and part of the filter-input computation 
section are bypassed. Otherwise, a test is made on INGUID to determine which set 
of guidance variables have been supplied (see Table 2-3). If |lNGLIDl=2, the roll 
angle about the body x-axis is obtained from the environment, provided INSENS is 
not zero. If |INGUID|*1. the roll angle Vo out the air-relative velocity vector is 
computed using air- relative velocity, position, and the direction cosine matrix. 

The command processing section simply sets the guidance- commanded pitch- 
rate (cjq) equal to zero when |INGUID|-1 and sets the guidance -commanded anglo- 
of-attack (oq) equal to the current measurement (a) when | INGUID} *2. 

The automatic mode filter inputs (errors in the angles of attack and sideslip, 
the appropriate roll angle, and the body rates) are then computed. Altitude and 
velocity magnitude filter inputs are either computed or zeroed, based on the value 
of MODE. These variables are used for automatic speed trim and aided guide -slope 
following during a manual approach. The computation of the desired values of 
altitude and velocity magnitude, H c and Uq. is described in Section 3C. 

When the Terminal Guidance Phase starts. INGUID changes from 1 to 2 and 
the set of guidance commands which is to be processed by the DFCS changes. When 
this occurs, the quantities * p , q i9Uin . and q osum used in the cruise filter update 
are set (initialized) to 0. 

Finally, a check is made to ensure that INGUID is negative. This avoids re- 
processing non-updated G&N inputs in subsequent DFCS passes. Control is then 
returned to the DFCS executive routine. 

The relatively simple Input Interface Routine presented here will undergo 
substantial change for the reasons which follow. The introduction of sensor models 
will directly affect the data input section of the interface. Future developments in 
the DFCS filters may dictate coordinate transformations in the computation of 
filter inputs. Finally, the simple command-processing section will employ buffer 
cells to protect the DFCS observations of guidance commands from asynchronous 
alteration. 
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B. Air Data Routine 


The Air Data Routine is called by the executive immediately after the Input 
Interface Routine during the atmospheric flight phases. It computes several para- 
meters required for aerodynamic control* 

The inputs to this routine are the following measurements which are computed 
by, or transmitted through, the Input Interface Routine: 

T A — dimensionless temperature 

r\ 

a — angle-of- attack (radians) 

H - barometric altitude (ft) 


U 

«dp 

The outputs are: 


air-relative velocity (ft /sec) 

2 

dynamic pressure (n/m ) 


mean (trim) angle-of- attack (radians) 


oi no •» , 

^prc^mations to the sine (cosine) of 

COSQ M 

U„, - mean airspeed (ft/sec) 

M 


Jfi - Mach number 

2 

q d p - dynamic pressure (lb/ft ) 

The means (which are the outputs of first-order filters with time constants 
equal to the slow sampling period) and the sine and cosine are updated every pass. 
The means are set equal to the values of their respective measured inputs in an 

initialization pass. 


The other outputs are only computed on slow -sample -rate passes, 
putations are shown in Fig. 2-2. 


The com- 
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START 



Yes 


7/i = 


u 


M 


Inputs: T ^ , a 9 H. 
Saved from past: o 


g a*V t s 


1051 [HT A (1-{T A )] 
dp 


q dp * °- 0209 


c 


END 


Fig. 2-2 Computations in Air Data Routine 








C, Outputs to Actuators 

The DFCS issues commands to two sets of control effectors: the Attitude 
Control Propulsion System (ACPS) jets and the Aerodynamic Control System (ACS) 

surfaces. 

In the entry and transition phases, pure torque commands are currently 
issued to the ACPS about any or all of the three body-axes. A clockwise-positive 
convention is assumed. A capability to interface with a more realistic jet model 
is in development; the design for this jet-selection capability is described in the 

next subsection. 

Commands are issued to the ACS during the atmospheric flight phases - 
currently entry, transition, and cruise. In entry, commanded aileron position 
<fca r ), commanded elevator position (4 e c ), and commanded speed-brake position 
(fcsb r ) are computed (in radians) and transmitted. In transition and cruise, the 
additional output commanded rudder position (6r c > is transmitted. The assumed 
convention for 6a c is such that a positive aileron deflection produces a positive 
torque about the body x-axis; the magnitude of 6a c is equal to the commanded 
left-elevon defection minus the commanded right-elevon deflection, where posi 
tive eleven deflections are defined as downward. The assumed convention for 6r Q 
is such that a positive rudder deflection produces a negative torque about the body 
z-axis. The assumed convention for 6e c is such that a positive elevator deflection 
produces a negative torque about t ie body y-axis; the magnitude of ie Q is equal to 
one half the sum of the commanded left-elevon deflection and the commanded right- 
elevon deflection, where positive elevon deflections are defined downward. The 
assumed convention for 6sb c is such that a positive speed-brake deflection increases 
drag; the magnitude of 6sb c is equal to the angle between the two sides of the split 

rudder. 

All the aerodynamic surface commands are limited within the DFCS prior to 
being issued. The limits magnitudes are shown in Table 2-4. The aileron and 
elevator commands are treated separately even though they apply to the same pair 
of actuators, the elevons, and are, therefore, interrelated; this question will be 
dealt with in the future. 

Elevon commands are related to aileron and elevator commands by the 
following equations: 


6L C = 6e c + 

1 

2 


6R C = 6e C ' 

X 

2 

6a c 
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Table 2-4: DFC5 Limits for ACS Commands 


Surface 

aileron (differential elevon) 
elevator (tandem elevon) 
rudder 


Command Symbol 

6a c 

6 Cri 


6r C 


6sb, 


Limits 


Radians 

Degrees 

±0. 3492 

±20.0 

+0. 2619 

+ 15.0 

-0. 7845 

-45. 0 

±0. 3492 

+20.0 

+ 1.5 

+ 85. 9 

0 

0 


speed brake 


D. 


AC PS Jet Selection Subroutine 


This section describes a scheme for converting angular-rate-change and_ 
translation requests into timed jet firing policies for the ACPS (Attitude Control 
Propulsion System). This routine has been tested with the DFCS program. 

1) Thruster Geometry 

The thruster configuration currently incorporated in the jet selection routine 
is the North American Rockwell ATP design. 

As shown in Fig. 2-3, the thrusters for the ATP configuration are grouped 
in 16 locations. There are 40 thmsters in all, and these are distributed among 
the 16 locations as shown in Tables 2-5 and 2-6. There are 16 thrusters in the 
nose of the fuselage and 24 in the "tail". 

2) Jet Selection Logic 

The DFCS Control Routine determines a commanded velocity change and calls 
the jet selection subroutine, which handles angular velocity components and trans- 
lational velocity components simultaneously with a single algorithm. A flag can be 
set instructing the jet selection subroutine to ignore translation, simplifying and 
speeding up the jet selection procedure. Furthermore, the jet selection logic can 
be instructed to deal with jet clusters or with the individual thrusters. 

Figure 2-4 is a functional flow diagram of the jet selection logic. A dis- 
cuss ion of this logic follov/s. 

The ACPS jet control authority subroutine is only fully executed when the 
jet selection routine is called for the first time in a mission phase. The control 
authority subroutine calculates the ACPS control authorities on the basis of total 
vehicle mass, the vehicle inertia matrix, and ^e center of mass location. The 
control authorities are calculated either for the jet clusters or for the individual 
thrusters, as is appropriate. The control authority is the angular and transla- 
tional acceleration in vehicle coordinates generated by firing a jet or cluster of 
jets. In addition, a cost per unit firing time is assigned to each jet. Currently, 
ail these costs are equal. When a jet is disabled (or restored) and the jet- 
selection logic is operating on clusters, the jet control authority subroutine is 
partially executed, appropriately altering the corresponding acceleration vector 
and cost element. 

The jet selection subroutine next checks its inputs to see whether or not 
there is a new velocity change request. If not, the jet selection subroutine call 
must be for the purpose of turning off jets which have completed their firing 
intervals. If the velocity change request is a new one. the linea" programming 
subroutine is called, which expresses the velocity change vector as a linear 
combination of jet acceleration vectors, the linear coefficients being the firing times. 
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Table 2-5: Distribution of ACPS Thrusters 


NR ATP Baseline Orbiter 


Thruster Location 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


dumber of Thrusters 

Thruster ID 

2 

10, 12 

2 

14, 16 

2 

13, 15 

2 

9. 11 

2 

6. 3 

2 

5, 7 

4 

17. 19, 21. 

4 

18, 20, 22, 

3 

32, 34, 36 

3 

26, 28, 30 

3 

25. 27, 29 

3 

31, 33, 35 

2 

2, 4 

2 

1. 3 

2 

38, 40 

2 

37, 39 


40 total 
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Table 2-B Pararie^ers of Individual fets 


ID Jet Locatio n - Body Station (inches) L h r £m£o nentg 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 
?0 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


318. 0 

- 96.0 

374. 0 

318. 0 

96. 0 

374. 0 

318. 0 

- 96. 0 

349. 0 

3 Id. 0 

96. 0 

349. 0 

318. 0 

- 109. 0 

368. 0 

318.0 

109. 0 

363. 0 

318.0 

-109.0 

355.0 

318.0 

109.0 

355.0 

3 id. 0 

-109.0 

346. 0 

318.0 

109.0 

346. 0 

330.0 

-109.0 

346. 0 

330.0 

109.0 

346.0 

318.0 

- 80.0 

367.0 

318. 0 

80.0 

367. 0 

330.0 

- 80.0 

367.0 

330.0 

80.0 

367. 0 

1538.5 

-158. 0 

425. 0 

1536.5 

158.0 

425. 0 

1532. 5 

-158.0 

438. 0 

1532. 5 

158.0 

438. 0 

1528. 5 

-158 0 

450.0 

1528.5 

158. 9 

450. 0 

1525.5 

-158. 0 

464. 0 

1525. 5 

158.0 

464. 0 

1541. 5 

-130. 0 

480.0 

1541. 5 

130. 0 

480.0 

1555. 5 

-130.0 

484.0 

1555. 5 

130.0 

484. 0 

1567. 5 

-130. 0 

487. 0 

1567. 5 

130 0 

487. 0 

1546. 5 

-130. 0 

416. 0 

1546.5 

130.0 

416.0 

1559.5 

-130.0 

420.0 

1559.5 

130.0 

420.0 

1572. 5 

-130. 0 

428.0 

1572. 5 

130. 0 

428.0 

1579. 5 

-149. 9 

468.0 

1579.5 

149.0 

468.0 

1579.5 

-149.0 

455.0 

1579.5 

149.0 

455.0 


-1.0 

0. 0 

0. 0 

-1.0 

0. U 

0. 0 

-1.0 

0.0 

0.0 

-1.0 

0. o 

0. 0 

0.0 

1.0 

0.0 

0.0 

-1.0 

0. 0 

0.0 
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Fig. 2-4 ACPS Jet Selection Functional Flow Diagram 


2-18 





Generally, many different jet combinations could satisfy a Riven velocity change 
request; if there is more than one such combination, the total cost is minimized 
in the selection process. Any jet which is flagged because of failure or mission- 
related constraints is omitted from consideration. 

The selected jets are turned on (if they are not already firing), and all others 
are turned off and removed from the internal jet-off list. The interrupt clock is 
set to return control to the jet selection subroutine when the next jet- off is due. 
and the "new velocity request flag is zeroed. 

The interface parameters for the jet selection package are: 

INPUT: 

a) Velocity change command vector. 

b) Flag for new velocity vector. 

c) ACPS jet exclusion list. 

d) Translation inclusion flag. 

e) Phase flag. 

OUTPUT: 

a) Jet- on^off commands. 

b) Interrupt clock settings. 

c) Zeroed flag for new request. 

Figure 2-5a presents features of the linear programming subroutine, CCWLJS, 
and Fig. 2- 5b describes its logic flew. 
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Linear Programming Subroutine 


PI’ R FOSE: Represents input vector W as a Linear combination of Unearl> 

- independent columns from matrix B, choosing columns from B 

to minimize cost of solution. 

INPUTS: LROW - Number of components m V. (Dimension of requests) 

LCOL * Number of columns in B matrix. (Number of jets or clusters) 

W * Vector (size I .ROW) to be analyzed as a linear combination oi 

columns from B. 

COST - Array (size LCOL) of strictly positive cost coefficients, one 
associated with each B column. 

B - Matrix (LROW by LCOI.). Each column represents the 

acceleration associated with a one-second firing of the corre- 
sponding jet in vehicle coordinates. Components may include 
rotational or translational acceleration elements, or both. 

LALARM * Non-zero only if B. LROW, or LCOL have been changed since 
previous call. 

JFAIL - Array (size L.COL) of flags, I th element non-zero to indicate 
jet failed. 

PARAMETERS - SAVING, MAXPAS, MAXITR. TOLRNC. can be chosen to 

match precision of the compilation and accuracy of input data. 


OUTPUTS: BI - Transformation matrix (LROW by LROW) going from vehicle 

~~ coordinate to coordinates based on selected columns of B. 

V - Solution vector (size LROW) of coefficients for selected columns 
of B matrix, satisfying: BI W = V. 

CF - Vector (size LROW) of cost coefficients selected from COST 
array, corresponding to chosen columns from B matrix. 

Z - Cost of solution, satisfying: Z = CF • V, 

LQLIST - Array (size LROW) giving the indices and sequence of selected 
B columns. 

MCOUNT - Number of basis element substitutions performed in current 
execution. 


LALARM - Zero, unless one of two cases occurs: 

1) MCOUNT has reached its maximum. 

2) Some B column has a negative cost component in the COST 
array. 


Fig. 2-5a Features of Linear Programming Subroutine, CCWLJS 
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Linear Programming Subroutine (Cont) 


ASSUMPTIONS: 


1) Cost array entries are all strictly positive, i. e. , greater 
than zero. 

2) Negative LQLIST outputs and negative V components are to 
be excluded from any implementation of the solution. 

3) COST array entries are -small compared to 1000. 

4) BI matrix is never altered by any other program. 

5) LALARM will be input as non-zero whenever changes occur 
in LROW, LCOL, or B. 


ENTRY : Call CCWLJS (B. W> 

EXIT: RETURN to caller. 


Fig. 2- 5a 


Features of Linear Programming Subroutine, CCWLJS (Cont) 
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Common initialization. 



If cumulative round-off in BI 
is large or if external changes 

are flagged in B, LROW or LCOL, 
then do complete initialization: 

* 

Null pass counter. BI is zero, 
with main diagonal equal signum 
of^W (+ 1 for W(I) a 0, - 1 for WU) 
negative). LQ negatives show 
non- jets. High cost insures re- 
placement of non- jet basis elements. 
Positive LRCORD indices corre- 
spond to available jets. 


If current U7 is in th e s ame 
sector as preceding W, use 

old Bf . For components in 
which sign is reversed, negate 
corresponding basis element, 
negate signum record component, 
set processing flag. 


Fig. 2-5b Flowchart of Linear Programming Subroutine, CCWLJS 
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M counts number of passes 
on this execution. 


0 - M 


00 



Latest version of solution, in 
terms of current basis. 


Bypass further processing, if 

BI is healthy and W sector is 
unchanged from previous entry. 


If no solution is found after 
MAXITR passes, set alarm, 
exit. 


Y(I) is cost of generating I 
original (vehicle) basis vector 
as linear combination of current 
basis vectors. 


Initialize the search for next 
jet to be included in basis. 
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Flowchart of Linear Programming Subroutine, CCWLJS (Cont) 






: •*»•.«**♦*» *W» m** 


I = 1 to LCOL, then go to 300 


Search for next jet to be 
introduced to basis, bypassing 
any which are failed or already 
in basis. 


LRCORD(I) <r 

or 

jet failed^ 


ves Bypass this 
jet 


GO TO 
200 


X » Y • B( I) - COST(I) 


X is the price penalty for sim- 
ulating a one- sec. firing of I th 
jet, using the basis rather 
than the jet itself. 


X > H 


X — H 
I - K 


Greatest saving per second is 
H; corresponding iet is K. 


GOTO 

L 200/ 



GO TO 
450 . 


If no jet was selected for 
inclusion, then exit. 


K m iet is to be included. Y 
is K* h jet acceleration vector 
described in current basis. 


Fig. 2- 5b Flowchart of Linear Programming Subroutine, CCWLJS (Cont) 













r I = 1 to LROW, then go to 330 




Choose T h member for replace 
ment so that new V( 0 is s 0, 
and no other V component 
changes sign. 


If no basis element can be 
legitimately replaced, set 
alarm, exit. 


Fig. 2-5b Flowchart of Linear Programming Subroutine, CCWLJS (Cont) 







Bad result. 


Fig. 2- 5b 


Flowchart of Linear Programming Subroutine, CCWLJS (Cont) 







EXIT 



For 1 = 1 to LRC'W then go to 460 
— — r 



Calculate cost of solution: 

Z = CF • V 


Negative firing times represent 
requests which cannot be met 
with currently non-faiied jets. 



Copy list of selected jets, 
alarm, and iteration count for 
external use: save internally. 


* 

NPASS reflects quality of BI 
matrix. 


Return to caller. 


Fig. 2-5b Flowchart of Linear Programming Subroutine, CCWLJS (Cont) 
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E. 


AC PS Jet Selection Logic for the ATP Conf iguration 


A jet selection routine is presented in this section that is more restricted 
than the design of the preceding section. It is essentially a table look-up; 
the firing policy for each possible command, as well as the individual jets 
an 4 the control authority corresponding to each jet cluster are pro- determined 
and stored for access by this on-line routine. The North American Rockwell 
ATP design for the AC PS , which has been presented in Fig. 2-3 and 1 ables 
2-5 and 2-6, was assumed in the off-line computations. 

The flowchart of Fig. 2- 6a describes in detail the jet selection logic 
for the ATP configuration. Figure 2- 6b defines the flowchart symbols. A 
discussion of some of the features of the logic follows. 

General - To provide attitude control with no translational disturbances, jet 
couples must be used. However, pure jet couples exist in only one axis, the 
roll axis. Jet couples to be used for pitch control (e. g. 3 & 12) and yaw 
control (e. g, 5 & 7), also effect roll. To achieve pure pitch or yaw two 
thruster couples must be fired to null the roll disturbance; for example, 3 and 1 
12 with 2 and 10 for pitch and 5 and 7 with 3 and 1 for yaw. Hence, couple 
control can require as many as four thrusters to fire whereas non- couple 
control may require only one thruster firing (e. g. yaw control). Additionally, 
the multiple thruster firings required for couple control produces large minimum 
impulse rate changes and these cause fast limit cycles and inefficient propellent 
usage. This logic has the option of attitude control with either coupled or 
non- coupled jets. 

Nose versus Tail Thrusters - for Non- coupled Jet Operation - The logic allows 
the crew to select either nose or tail thrusters for control of each axis. 

This serves two purposes 

a. Comparing nose thrusters to tail thrusters, the nose thrusters have 
shorter roll moment arms, but longer pitch and yaw moment arms. Hence, 
the crew can select the preferable moment arm length for each axis. 

b. This is a simple means of accounting for thruster iailures. Should 
all the thrusters in any single location fail, control can still be maintained 
(Note: Individual failures in any jet location are accounted for by placing another 
jet of the same location in "common" with the location, see page 16 of the 
flowcharts). 

Phase Plane Jet Selection Logic Interaction - (Non- coupled Operation) - 
Thrusters in locations 1 through 4, and 9 through 12, produce both Ditch and 
roll torques. Hence, it takes two thrusters to provide either pure roll or pure 
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pitch. However, smaller minimum impulse rate changes are possible by firing 
only one thruster. Therefore, even if there is a phase plane command from 
only one axis (roll or pitch), the logic will attempt to do the le-er (fire only 
one thruster) if the resulting disturbance in the "non-command axis reduces 
the rate error (called RRy or in ^ ie 


Use of the tail thrusters for yaw control causes a roll disturbance: should 
the roll disturbance exceed a rate limit, the logic will consider the disturbance 
part of the phase plane command and cause corrective roll thrusters to fire. 


Translations - No combination of thrusters will produce V or Z translations 
with negligible rotational disturbances. Thrusters must be modulated to null 
rotational disturbances. This modulation could be performed either open 
loop in the jet selection logic, or closed loop through the phase plane logic 
(the crossing of a switching line) and jet selection logic. The latter method 
was chosen because it will not produce unnecessary modulation (thruster 
firings) during small Y or Z translations (i. e. . when rotational disturbance 
is not long enough to cause switching line crossing). 


During +X- translations, it is poss u e to control yaw by modulating the +X 
thrusters; this option is available i n the logic. 


Jet Coupling and Logic Flow - With respect to thruster control actions, there 
are two basic types of thrusters 

a. Those that produce Y translation and yaw. 

b. Those that produce Z translation and roll and pitch. 

Sequentially, the logic first calculates commands for Y translation and yaw and 
then issues commands for Z translation, roll, and pitch control. 
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Fig. 2-6 b - DEFINITION OF SYMBOLS FOR JET 
SELECTION LOGIC 

Ccnnanda for translation, either positive, none, or 
negative, depending on +1, 0, or -1 

C csnanded rate changes frcm switching logic. 

mi 

Angular Acceleration for i jet for either roll, 
pitch, or yaw. 

Value of 1 (not 0) signifies request for long maeent 
an thrusters 

Indices for jets 

Min lmm impulse time 

Roll rate disturbance for yaw end/or Y translation firings 
Roll rate limit before HDFY is added directly to 0 C 
Roll rate error from phase plane logic when 0., is aero 
Index indicating presence (l) or absence (0) of roll caranjand 
Sample time 
Jet Indices 

Maximus roll compensation time available by unbalancing 
pitch firing times 

Snail master to account for round off in ccaiperisions 
(say 10"**) 

One jet roll correction time 
Ret roll comsend, i.e., HR ♦ 0 C 

# 

Roll rate limit before RR^is added to 0 C during translation 
Sample time plua E 

4 

Pitch rate error from phase plane logic when ©^ is aero 
Couple option 
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SECTION 3 


DECS - CRUISE PHASE 
by 

Albert G. Engol, Franklin H. Moss, Alexander Penchuk 


A , Introduction 

The crui 3 ing flight phase of an orbital Space Shuttle flight begins at Mach 0.9 
and ends at touchdown; it follows the high-angle -of-attack.- high-Mach-number 
entry and transition phases. In the near future, the cruise phase will be extended 
to deal with Shuttle powered horizontal flights from lift-off to touchdown. During 
this phase the control is accomplished by moving the aerodynamic control surfaces: 
elevons, rudder, and rudder flare (speed brake). 

For the cruise phase the DFCS contains four modes of operation: three of 

them are manual modes and one is automatic. The manual modes are Direct, SAS 
(stability augmentation system), and RCAH (rate command attitude hold). The 
automatic mode is subdivided into two regions of operation which are determined 
by the guidance. The first of these using "approach ' guidance extends from an 
altitude of about 40000 feet down to roughly 7000 feet. Terminal guidance is used 
from 7000 feet to touchdown. Appropriate control laws are used to accommodate 
these guidance modes. 

Stabilization of the vehicle in the cruise pnase is done with digital compen- 
sation. The compensation filters are provided with two sets o; gains, a fixed set 
and a variable set. The variable set is used to account for variations in flight 
conditions and, consequently, in vehicle dynamics. The gains are changed as a 
function of dynamic pressure. Ali gains are based on a sampling rate of 10 
samples per second. 

The overall DFCS structure and the interfaces have been described in 
the preceeding two sections. The description of the cruise phase auto- 
pilot is presented in this section in the following sequence. First, the essential 
features of the cruise phase control modes are described with simplified functional 
block diagrams . These diagrams show all the feedbacks and compensations used 
at this time but exclude functions which are provided in the program but are 
effectively disabled by gain settings. The aided- manual- approach feature is 
also described. The digital filters are then discussed, followed by a presentation 
of all the control gains and filter gains, including a table of fixed gains and graphs 
of those that are varied with dynamic pressure. Finally, the actual detailed 
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implementation in the DFCS is presented in two forms: detailed block diagrams 
of all the control modes, and complete flowcharts of each of the major subroutines 
Recently developed design improvements are presented in an appendix. 

Additional information about the cruise phase autopilot can be found in 
References 3-1 through 3-3. 

B. Control Modes 
1) Direct Mode 

in the direct mode, pilot s commands are passed by the Dfr CS directly to the 
surfaces, but with a provision for crossfeeding a scaled speed brake command to 
the elevator to offset the resulting pitching moment. Also, a crossfeed from 
aileron to rudder is provided to offset a possible adverse yawing moment due to 
aileron deflection. 


r 



ELEVATOR 

COMMAND 


SPEED BRAKE 
COMMAND 


MANUAL 



RUDDER 

COMMAND 


AILERON 

COMMAND 


Fig. 3-1 Simplified Functional Block Diagram of the Direct Mode 
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2) SAS Mode 


In the Stability Augmentation System mode two rate feedback loops with 
simple gain compensation are employed. These are the pitch-rate and roll-rate 
loops. A pitch rate bias proportional to the yaw rate and the (estimated) roll angle 
is added to the pilot command to prevent adverse nose -down response due to steady 
turn rates and a crossfeed from the speed brake to the elevator is also provided. 
The commands to the rudder and speed brake go unchanged to the servos. 


MAN! A I. SPEED 
BRAKE COMMAND 


M \N! \E 
PITCH HATE 
t f >MM AND 


MAN! XL 

rc » u. k\ n: 

Ct >\! M AND 


MW \E 

H* ODER ( < Pj MAN 



<PEEI» BRAKE 
COMM AND 


V\ K\ \ Ji *R H »M\I \\L> 


m i; vs i k i ; i > 

prrt »ih \ n: 


\? \ I'Mi »\ < ( » J M \ Mi 

ME \Sf HE!) 

H( HA. H vn: 


WEXSI RED 
‘i U\ H \ I E 

H I DDE!? < < ‘MM \NI) 


Fig. 3-2 Simplified Functional Block Diagram of the SAS Mode 


?) HCAH Mode 

The Rate -Command- Attitude- Hold mode provides two rate control systems 
(pitch rate and roll rate), each of which has proportional-plus-integral compensa- 
tion. In addition, yaw rate, filtered by a "washout 11 filter, and sideslip are used 
as feedback quantities to the rudder. A pitch rate bias and a speed brake cross - 
feed are added into the elevator channel, as they are in the SAS mode. The man- 
ual command to tne speed brake is passed unchanged by the DFCS. 
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Fig. 3-3 Simplified Functional Block 
Diagram of the RCAH Mode 


Speed 

BraK*' 

Command 


Elevator 

Command 


Measured 
Pitch Rate 


Aileron 

Command 


Measured 
Roll Rate 


Rudder 

Command 


Measured 
Yaw Rate 


Measured 

Sideslip 
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4) Automatic Mode 


The automatic mode is subdivided into two regions of '*pera 
the following table: 


ion according to 


Name of 
guidance phase 

Apn roach 
I er mi na 1 


Approximate Guidance 

altitude range commands 


40 nno i.o Toon feet 



7 i *00 feet to touchdown 


q G . 


0 


b 


G 


During the approach phase the longitudinal control is of the angle of attack. In the 
terminal phase, longitudinal control is accomplished by the pitch rate loop which 
is identical to the RCAH loop described previously (including pitch rate compensa- 
tion during steady turns). In both cases the lateral control is of the roll attitude. 
However, the roll angle commands are different in the two guidance cases: roll 
angle is taken about the velocity vector in approach and about the body (vehicle 
x-axis ) in terminal phase. For both cases the form of the lateral control system 
and the gains are the same. 

The rudder loop is implemented similar to the RCAH mode except that a 
provision is made to accept yaw rate and sideslip guidance commands for automatic 
decrab maneuver. Rudder pedal command provides for this maneuver to be made 
manually. 

Inputs to the filters in the automatic modes are formed in the input interface 
subroutine by summing the guidance commands with appropriate quantities derived 
from measurements. 


C. Aided Landing Approach with ’’Canned Trajectory 

During landing approach, the DFCS can assist the pilot by automatically 
following velocity and altitude profiles which are stored functions of range-to- 
touchdown; hence an external source of range-to-go must be provided for this 
feature (which is switched on or off by a flagword). The control loops for the 
aided -landing-approach mode are the manual SAS and RCAH modes described above. 

Two trajectory curves are stored, velocity versus range and altitude versus 
range. The velocity trajectory is followed by use of the speed brake (or rudder 
flare). The altitude trajectory is followed by use of the elevator. Pitching moment 
caused by rudder flare deflection is countered by an additional deflection of the 
elev tor. The range-to-go is broken into four segments. Trajectory information 
is stored at the endpoints of each segment for velocity and altitude. The desired 
velocity and altitude at the endpoints of the segments is stored in the Pad Load (see 
Fig. 3-5). 
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Fig. 3-4. Simplified Functional Block Diagram of the 
Automatic Mode. 
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ra lector y Values 



Fig. *-5 Velocity and Altitude Profiles of the Canned Trajectory' 


Kor a ;,. tu , S ... th. segment-. a.-: Umlin* flan* (ending m touch, Urn n * Uh a 

.ink rate). straight line gl*'**’ • approach pulM-p, an,! -.t, > |n i K UtU si i * ^ 

-truight line. The flare and pullup ailments for iltitud.- are each iUHc ,unc 
H . such that the a lope ami position the curve cu.km!, with the adiacerU an.u t ht 
in*, segments, Coefficients fur th. *■ mte relation, are , ompu.ed in the OK s in 
the Initialization Routine. 


The three velocity segment.-, of the reference la 
touchdown are straight lines. The fourth segment u 
linear IAS (indicated airspeed) versus range. 


luting trajectory 
exponential and 


nearest 
results in a 


R . is the rango-to-go (to touchdown). When INSENS is set to unity, the DFCS 
computes V. % from data from the guidance and navigation (GAN) program. The DFCS 
determines H aegment the vehicle is in by comparing K GO with the five .tored .values 
of range. Having found the segment the program proceeds to compute U Q and H r> . the 
velocity command and altitude command. These commands are differenced with the 
estimated velocity and altitude in the Invut Interface Routine and the results are 
filtered. The filter outputs arc used in the Control Routine to determine the speed 
brake command and terms in the elevator command. 


An additional feature exists which adds a constant command bias to the speed 
brake and a corresponding bias to the elevator within each of the four trajectory aug- 
ments. The values of these biases are to l*- determined. 

The aided -landing commands are summed with the manual commands. In 
principal, the aided -glide -slope feature could (with higher gams) make an automatic 
landing without manual intervention; however, the intent is to use relatively low 
gams in the glide-slope feedback, resulting in long-term flightpath corrections 
which do not compel* with the pilot’s short-term commands. 


When R G( is greater than R 5> the fl»gword which switches aided- landing - 
approach on or off is used to enable the Mach-trim feature. Nominal values of 
elevator deflection are stored as a function of Mach number: these values 
approximately null the sum of aerodynamic pitching moments on the spacecraft, 
reducing the need for manual trimming as Mach number decreases. Tne trim 
value of elevator deflection is obtained by linear interpolation from a stored table: 
the values in this table are to be determined. 

D * Digita l Filters and their Implementation 

Control filters for the cruising flight mode are classical, recursive. 2 n -order 
digital filters, whose z transform is D(z) : 

-1 -2 

a l * a 2 z * V 
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The z** ..alculation blocks are achieved in the DFCS programming by the 
relative sequence of instructions, which permits ihe DFCS cycle time lag to occur 
at the point sc designated. In the coding, each state filter is divided into two parts 
in order to minimize the computational transport lag between data input and control 
command. The first part, which occurs in the Filter Update Routine <«=* ig. 3-16) 
is limited to updating the filter estimate according to the most current, measurement 


x . * a.x, + y 
out 1 in 1 gum 





The Filter Update Routine is sequentially followed by the Control Routine 

(Fitf, !-t7) which issues control commands based upon the most current filter 

output value x , Filter data ’pushdown” ts then accomplished in the Filter Push 
1 out 

down Routine ( Fig. 1 1 H) by the following two equations: 


x 

o 


sum 


a 4 x in * Vout 


This sequence of instructions requite in the following equivalent difference 
equation for filter output* 


out 




X. 

in-2 


■ 


out.} 5 out. 2 


where subscripts -1 and -2 denote past values, one and two samples ago. 

Control loops which utilise 2 nd -order filters are available for control of the 
following quantities: 

Manual Modes: roll rate 

pitch rate 
yaw rate 
sideslip angle 
normal acceleration 
roll rate to rudder 

Automatic Modes: roll angle 

pitch angle 
yaw angle 
angle of attack 
sideslip angle 

Manual and Automatic Modes: earth -relative velocity 

altitude 

All rate and angular control functions are performed at the fast sampling 
rate (10 sps). velocity control at the medium sampling rate (2 Bps), and altitude 
control at the slow sampling rate (1/2 sps). 

A recursive, 1®- order washout filter is used in the yaw rate feedback path 
for rudder control; * transform of this filter is: 


D(z) 


OjjO-i- 1 ) 


1 - G 37 ! 


-1 


3 -1C 


A block diagram of the update and pushdown equations for the washout filter are 
presented below. 



FU«.r update: r >>(h • 0,,*^ ♦ r >um 

FU«er puahdown: r lum • -G,,!-^ * G 37 ’ r w „ h 

A control bias for the maintanence of trim in manual -mode operation ia 
achieved through the integration of discrete signals generated by the trim buttons 
on the pilot’s hand controller. The z transform of the integrator is: 




A block diagram of the update and pushdown equations for the pitch channel are as 
follows: 



Filter update: 


**«r* ( K 15* T F“'.r)/ 2 
^tr ~ ie tr + 'krsum 


Filter pushdown; sum = « <* tr + 
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E. Gain* for the Cruise Phase 

A list of fixed gain# for the cruise phase manual and automatic modes is 
given in Table 3-2. The gains are re-evaluated every time a switching of modes 
takes place. These gains yield good performance in all modes of operation for 
flight conditions where dynamic pressure la about 225 Ibs/sq ft. 

Variations in flight conditions make it desirable to change certain gains contin- 
uously during the flight. This ir done in the Parameter Estimation Routine (Fig. 1-1*0 
every alow pass if the parameter estimation is enabled. Fourteen gains are changed 
as a function of dynamic preaaure, as shown in Fig. 3-9 and 3-10. The gains were 
computed and stored for seven dynamic pressures (100.150, 225, 300, 425, 525 and 
850 lbs /aq ft). A linear interpolation is made between these points. Outside of the 
dynamic pressure range, the gains stay equal to the boundary values. 


Table 3-2: Fixed Gains for Cruise Phase 
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Scheduled Gains See Fig. 5-9 and 5-10 

+ In SAS mode, these gains are: G., - G 2 + Ky ; K fl = O.o ; >• 0.0 ; Kj Q * 0.0 

$For alpha (approach) gu^ii-.vc. K q « 0.08*) 
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a) Automatic Mode, Alpha Loop (Approach Guidance) Gains 
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d) Pitch- Rate- Loop Gain in SAS Mode 


Fig. 3 -ft Cruise Phase Longitudinal Control Gains 
Versus Dynamic Pressure (Cont) 
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c) Gains Used in Automatic Mode Only 


V'w« 2 ' 


Fig. 3-10 Cruise Phase lateral Control Gams Versus 
Dynamic Pressure (Cont) 
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F. Detailed Block Diagrams of Cruise Contro l Modes 

This section presents .» detailed control-engineer's description <>f all control 
modes of the cruise phase by means of functional block diagrams. 1‘hese include 
all feedbacks, crossfeeds. limiter t. switches, as well as compensation filters in 
z-transform notation. The gains are properly subscripted and their values can be 
obtained from Table \-2 and Fig. i-9 and '<-10 on the proceeding pages. 
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-12A. Block Diagram of the Aileron SAS/HCAH Loops 
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Fig. 3-12B. Block Diagram of the Elevator SAS/RCAH Loop 































-13B. Block Diagram of the Elevator Auto Loop 



















G. Flowcharts of Cruise Phase Computations 

The flowcharts in Fig. 3-16 through 3-19 are a detailed presentation of the logic 
in cruise phase digital autopilot. They describe the Filter Update. Control. Filter 
Pushdown, and Canned Approach Trajectory and Parameter Estimation routines respec- 
tively. The numerical values of the control gains and filter gains can be determined 
from Table 3-2 and Fig. 3-9 and 3-10. 

Fig. 3-20 shows the special logic exercised in an initialization pass m the 

cruise phase. 
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Fig. 3-18 Filter Pushdown Routine for Cruise Phase 
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Fig. 3-20 Initialization Logic for Cruise Phase (Pad Load Manipulation Routine 
and Filter and Parameter Initialization Routine). 
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APPENDIX 3 A 


Preliminary Designs for Cruise Phase, Incorporating Powered 
Horizontal Flight Requirements 


This appendix documents a number of changes and additions to the cruise de- 
sign; many are for use particularly in powered flights (ferry flights and develop- 
mental powered horizontal flight tests). These designs are in a state of flux. 
However, they are presented at this time in order to point out areas that need 
attention and, in some cases, to present preliminary designs. 

The algorithms for control of powered flight are integrated into the cruise 
phase of the DFCS. It is recognized that decisions will have to be made with re- 
gard to partitioning when the software for the first horizontal powered flights is 
prepared, in order to ensure that the continuing development for the other mission 
phases does not impact the software that will be executed in the particular flight. 

Table 3A-1 presents an overview of the Powered Horizontal Flight (PHF) 
DFCS modes of operation. Manual operations imply active pilot control of all outer 
and some inner- loop variables, depending upon mode, and do not differ significantly 
from their unpowered counterparts. Flight Director operations also imply active 
pilot control of outer- loop quantities, but the pilot is aided by flight director fly- 
to" indications of desired control actions . Aided Glideslope operation provides 
low-gain following of a canned vertical- plane trajectory, but the pilot must fine- 
tune the trajectory with manual inputs; also, he must perform all lateral control. 
Semi automatic operation gives the pilot "hands off" capability with respect to 
selected outer- loop closures; he must configure the autopilot, and control man- 
ually those flight variables for which he has not (or cannot) select semi-automatic 
control. Fully-autormtic control implies "hands-off" operation. Commands to the 
DFCS are generated by sophisticated guidance /navigation modules. 

Manual operation takes place in one of three DFCS modes. In direct (DIR) 
mode (Fig. 3. A-2), inputs from the stick/pedal and from the speedbrake /throttle 
qradrants are sent directly to the appropriate control actuation system. Beep- 
trim" is available for stick and pedal inputs. A crossfeed from aileron command 
to rudder command provides a level of turn coordination consistent with this modo. 
Similar crossfeeds from thrust and speedbrake (manual) commands are used for 
elevator trim biasing (Fig. 3. A-14i. 
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In the stability augmentation modes tSAS, Figures 3.A-3, -4, -5), controlled 
variable (or damper Iood) feedbacks are introduced as follows: 

elevator loop - pitch rate 
aileron loop - roll rate 

rudder loop * yaw rate (after high-pass washout filtering) . 
lateral acceleration (or sideslip) 

The elevator and aileron loops are provided with forward- loop gain; thus the desig- 
nation "rate command" (RC) is appropriate for these loops. Crossfeeds are pro- 
vided as follows: 

elevator - thrust and speedbrake as in DIR, mach trim, yaw rate 

multiplied by an approximation of tan $ fcr turn-coordina- 
tion 

aileron - yaw rate and lateral acceleration (or sideslip) 
rudder - roll rate error for turn coordination 

Command augmentation is provided in the Rate Command Attitude Hold (RCAH) 
mode. Forward loop integral compensation assures unit-gain rate-tracking in the 
pitch and roll loops. The yaw loop retains its SAS configuration (unit- gain with no 
integration in the forward loop). 

C* command augmentation (CAS) is now provided. The controlled variable 
is a blend of normal acceleration and pitch rate: 



where 

g = gravitational acceleration constant 

The weighting factor V co must be empirically determined; it represents the airspeed 
at which the acceleration and pitch rate components contribute "equally" to the C 
feedback variable. The C* response to pitch stick deflection is relatively invariant 
over a wide range of flight conditions, compared with pitch rate alone, and pro- 
vides a natural response to pilot inputs, being sensitive to pitch rate at low vel- 
ocities and to normal acceleration at high velocities. 

Manual modes (DIR. SAS/RC. RCAH) are intended for backup operation where 
outer-loop controlled- variable measurements (♦, 6, <1, H, U, etc.) are not readily 
available. Loop configuration and performance are intended to be typical of that 
found in modern transport aircraft. The F-8 Digital-Fly-By-Wire Program 
(NASA/FRC) has already demonstrated digital fly-by-wire flight control using 
these manual modes of operation, and the current design to a large extent reflects 
findings from that program. 
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Flight Director (FD) operation provides numerous outer-loop capabilities, 
but the loops themselves are closed by the pilot, flying in one of the manual modes 
(DIR. SAS/RC. RCAH). The signals which drive the FD displays/ needles are all 
available as error signals n the appropriate semi-automatic loops. Several 
options are being considered, as follows: 

pitch bar - pitch altitude error (» in > 
roll bar - roll attitude error (0 in > 
high/low needle - altitude error (H in > 

altitude rate error (H in > 

1LS glideslope error (e^g) 
left/ right needle - lateral displacement error (Y. n > 
heading error 
VOR error < e V Qj^ 

ILS localizer error ,e L oC* 
fast /slow indicator - velocity error <U in ) 

All such displays will indicate "fly-to” directions of desired response. 

Sensitivity of the displays will be individually controllable. 

Aided Glideslope (AGS) operation (Fig. 3. A-19) provides low-gain altitude 
loop closure through the RC/RCAH pitch rate loop, and low -gain velocity loop 
closure through either the automatic speed brake loop (AUTOSB, Fig. 3. A-ll) or 
the automatic throttle loop (AUTOTH, Fig. 3.A-13). Details of the AGS 
trajectory are given in Section 3.C. Since the loop closures are low-gain, the 
pilot must supply whatever additional stick/ pedal and speedbrake/ throttle over- 
ride commands are necessary to trim up the trajectory. In particular he must 
provide all lateral axis control, 

Semi-Automatic operation provides full autopilot outer-loop closure capabilities. 
Several options are provided, permitting individual selection or reasonable com- 
binations thereof. Loop gains are high, and with respect to the option(s) chosen the 
flight is "hands-off". Variables left uncontrolled are the responsibility of the 
pilot, using RC/RCAH modes. Pilot override is always available. Disengaging of 
outer-loop features, as in control stick steering, has not yet been addressed. 
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Semi-Automatic options fall into several categories, as follows: 

ATTITUDE CONTROL (*, *. outer loops ) 

ALTITUDE CONTROL (H, H outer loops) 

SPEED CONTROL (AUTSB or AUTOTH loop closure 
around indicated air speed, inertial speed, 
or ground track speed; pitch rate loop closure 
around Mach) 

RADIO AIDS Tracking (VOR and ILS GS/LOC capture /track) 
3D-GLIDESLOPE (canned trajectories) 

ALIGN /Decrab (<4, r, Y loops) 

RNAV (Simple trajectory- following via way-point segments) 

Selectable options under these categories are given in Table 3. A-2. Select 
(SEL) options imply pilot-insertion of data, whether via thumbwheel (©increments 
or * increments), dedicated window, or general purpose keyboard (altitude or 
altitude rate desired, etc. ) The options indicated herein are commonly found 

on commercial and military aircraft (B747, L1011, DC10, B70, F8. etc.)- 

Fully-automatic "hands-off" operation is expected for the DFCS automatic 
(AUTO) mode (Figures 3.A-6 through 3. A-8). A number of interfaces are 
provided in each axis, including direct guidance inputs to the elevator, aileron, 
rudder, and speedbrake or throttle commands. Primary guidance interfaces are 
as follows: 

Speed: velocity (U^.) 

Pitch: a- guidance (Oq) or Q guidance (q^) 

Roll: f v ~ guidance or ^-guidance 

Altitude: altitude (H Q ) 

The AUTO mode loops remain relatively unaffected by the addition of semi- 
automatic and flight director capability; detailed description is given in the main 
body of Section 3. Speed and altitude loop closures are handled differently how- 
ever; see Figures 3.A-9, 3.A-10, 3.A-12, and 3.A-17. 

In summary, the baseline DFCS design forreturn-from-orbit cruise phase 
remains relatively unchanged. Significant additions and modifications shown or 
implied by figures 3.A-1 through 3.A-20 are as follows: 

pilot input processing: bias removal, deadband de sensitization, 
and prefiltering (Fig. 3.Arl) 
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generalization of semi-automatic elevator auto-trim (Fig. 3. A- 14) 
speed control loop configuration (Fig. 3.A-9, 3.A-10. and 3.A-12) 
throttle and speedbrake control and interconnect logic, including 
manual and automatic throttle (MANTH and AUTOTH) 
logic (Fig. 3.A-11, and3.A-13) 
altitude control loop configuration (Fig. 3. A- 17) 

C* command augmentation mode (Fig. 3.A-20) 

semi-automatic and flight director operational modes (Table 3. A-2) 



N 


TABLE 3. A- 1 

Powered Horizontal Flight 
Digital Flight Control System 
Modes of Operation 


MANUAL OPERATION 
Direct (DIR) 

Stability Augmentation (SAS'RC) 

Command Augmentation (RCAH, C ) 

Manual Speedbrake (MANSB); Manual Throttle (MANTH) 

FLIGHT DIRECTOR OPERATION (manual modes) 

Pilot closes outer loops 

Error signals from semi-automatic loops displayed 

AIDED GLIDESLOPE OPERATION (SAS/ RC or RCAH) 

Low gain H- loop and U-loop closures 

Pilot responsible for fine tuning final trajectory 

Automatic speed brake (AUTOSB) or Automatic Throttle (AUTOTH) 

SEMI-AUTOMATIC OPERATION (SAS/RC or RCAH) 

Attitude Control Speed Control 

Altitude Control Radio Aids (VOR/ILS) 

3-D Glideslope Align (Decrab) 

RNAV Trajectory Following 

FULLY AUTOMATIC OPERATION (AUTO Mode) 

q-GUID or o-GUID; GUID or » x GUID 

MISCELLANEOUS OPTIONAL FEATURES 

Semi-automatic elevator trim 

Turn coordination (pitch rate and/or rudder) 

AUTOSB to MANTH or AUTOTH to MANSB interconnect 
Control stick steering override of semi-automatic control 
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TABLE 3. A- 2 


Powered Horizontal Flight 
Digital Flight Control System 
Semi-Automatic Control Options 


attitude CONTROL 

e- ATT, SEL/HOLD 
0- ATT, SEL/HOLD 
HDG, SEL/HOLD 

ALTITUDE CONTROL 

ALT. SEL/HOLD 
ALTRATE 
FLARE (H=-KH) 

SPEED CONTROL 

MACH SEL/HOLD 

SPD SEL/HOLD (IAS, INS, GTS) 

RADIO AIDS 

VOR CAP/TRK 
ILS GS CAP/TRK 
ILS LOC CAP/TRK 

3-D GLIDESLOPE 

H. U, Y vs. RGO 

ALIGN FOR LANDING 

$, r, Y CONTROL 


RNAV 


Way Points 

Simple Straight Line Trajectories 
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Fig. 3. A- 3 Aileron Loops for SAS/RCAH Modes 







































































Fig. 3. A- 11 Speed Brake Control Logic 
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13 Throttle Control Logic 










14 Elevator Trim-Bias Circuit 












Interface Signals for Semi-Automatic Roll Rate Command Logic 








Automatic Pitch Rate Command for Altitu* 













Interface Signals for Semi-Automatic Pitch Rate Command L< 
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20 Elevator Loops for CSTAR Option 
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Fig. 3. A- 21 Semi-Automatic HNAV Trajectory Following Logic 




SECTION 4 


TRANSITION PHASE 
by 

Richard D. Goss 


A. Introduction 

The transition phase of the return from orbit folio., the entry phase and is 
defined to be that port of the trajeotory for which Mach nutr.be- lies betweer o.O and 

0 . 9 . 

Attitude control during the transition phase is maintained by a combination of 
aerodynamic control surfaces and control jets. Yaw and roll control are achieved 
through the use of aileron (differential elevon) deflection, rudder deflection and yaw 
and roll control jets. Pitch control is achieved through the use of elevator (tandem 
elevon) deflection and pitch control jets. The use of yaw control jets is required*, 
maintain stability during the early phase of the transition trajectory when the high- 
angle-of-atUck/high-Mach-number operating conditions result in the rudder be ng 
ineffective. The roll jets and pitch jets act to support the elevons to prevent the 
buildup of high vehicle attitude rates. 

The control law for the aerodynamic controllers is a linear control law which 
is designed to control perturbations about the desired state. The control gains are 
scheduled as a function of Mach-number and dynamic pressure. The gams have been 
computed Tor a range of Mach-munber/ dynamic-pressure operating conditions by 
attempting to match closed loop system response to an "ideal" or desired system 
response for those operating conditions. 

The vehicle attitude and attitude- rate inputs to the controllers are provided by 
a linear state estimator (or filter). The filter is partitioned into longitudinal and 
lateral-directional modes. The longitudinal and lateral states contain only the com- 
ponents required to describe the dominant modes of vehicle flight. These states 
will be described as the reduced "states”. Additional state components used for 
extrapolation of the reduced state are obtained from guidance and navigation estimates. 
These components are described as the "observed states . 

Although the state estimator equations are incorporated in the current DFCS, 
the equations have been implemented to date only l.i a trivia) form. That is. gains 
have been selected which pass inputs through the filter unmodified. The present 
version of the DFCS assur.es a vehicle-body frame-of- reference for both angle and 
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rate estimation. The control laws, however, are designed to U3e inputs which 
describe vehicle orientation in an earth- relative velocity coordinate system and 
vehicle angular rates ir. a body coordinate system. The required angle inputs to 
the control laws are presently being supplied directly by the Input interface Routine 
described in Section 4 -2 A. 

The state estimator is ir. a developmental stage at present and is expected to 
change substantially. A determination of sensor requirements must be made before 
proceeding with filter redesign. The filter design is strongly dependent on the 
number of measured quantities which are available throughout the transition phase. 

If aerodynamic sensors are available, procedures for the most effective mixing of 
aerodynamic and inertial measurements must be developed. In addition, a deter- 
mination must be mode of the most effective coordinate frame in which to process 
filter inputs as well as the most efficient way to transform filter data to the form 
required by the control laws. 

A simple block diagram which provides an overview of control in the transi- 
tion phase is shown in Fig. 4-1. It illustrates the interconnections between the 
various major routines during this phase. The next four subsections of this document 
describe the algorithms which are implemented in the Filter Update, Control, Filter 
Pushdown, and Parameter Estimation routines respectively. Each of these descrip- 
tions consists of a brief explanation followed by a detailed description from the point 
of view of how the algorithms appear in the DFCS computer program. The final 
transition phase subsection presents another view: engineering block diagrams for j! 

each of the control channels from the Input Interface Routine to the setting of the j 

actuator commands. ; 

j 

B. Filter Update 

1 # Longitudinal Estimator j 

a) Engineering Description 

The reduced state of the longitudinal estimator consists of two components — 
pitch attitude and pitch rate. The reduced state measurements are incorporated in 
the linear estimator using sequential corrections. 

The observed longitudinal state components are angle-of-attack and velocity. 

In the current version of the DFCS these components are assumed to be obtained 
from guidance -and- navigation estimates or sensors with no additional filtering 
applied. In subsequent versions of the DFCS, additional filtering will be applied 
to the angle-of-attack estimate to account for high frequency vehicle attitude dyna- 
mics. 
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Fig. 4- 1 Overview of Transition Phase Control 








The longitudinal estimator also includes a pitch- rate compensation filter* 
This filter estimates the steady-state pitch rate induced in a coordinated turn. 
Since it is not desired to feed this pitch- rate bias oack to the control laws, the 
measured pitch rate is modified by this bias. The modified pitch rate is then 
received as an input by the linear state estimator. 


b) Im plementation 
Definitions 


Reduced longitudinal 
state vector: 


•o' 6 


- q 


(Pitch attitude error) 
(Pitch rate error) 


(9 and q are the actual values of vehicle pitch attitude and pitch rate, 

and 9 ~ and are the desired values from guidance of vehicle pitch attitude and 
G G 

pitch rate.) 


Reduced state measurement: 


9 in* 9 G- 6 


"in * “G ' " 


(9 . and q, are inputs to the filter update routine.) 
in nn 


Computations 

Pitch- Rate Compensation: 


= 0.999 d_ + T^ p 


q in ~ q in 


P F 
+ r ( d 


- *3 


P 3 P 


Notes: 


1 ) 


The term r (4 +4- ♦ 3 ) is the estimate of the pitch rate induced in 

pJp 

a coordinated turn 


2) (d + i *3 is an approximation to tan ( d ) 

pap r 

3) The factor T F is the fast sampling period (0. 1 sec) 

4) The factor 0.999 (which yields a time constant of 100) is a constant 
value in the current DFCS. 
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Sequential State Vector Correction (at Fast Sampling Rate) 
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where: 
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v out 
^out 
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and 


is the state estimate propagated from the previous 
sampling time to the present sampling time. 


is the state estimate folic wing the first correction. 


is 


the state estimate following the second correction. 


are the gain vectors which are multiplied by the measure- 
ment residuals. 


Angle-of-Attack and Velocity Updates 


a . * oi. 
out m 


u o„, * u ln 


(no filtering applied to a and U. a to be filtered at fast 
rate in subsequent versions of DFCS.) 


where: 


"in S *G “ a 
«in * U C ' 5 


2. 1 .atpral- Directiona l Estimator 

a) Engineering Description 

The reduced state of the lateral- directional estimator consists of four com- 
ponents-yaw attitude, yaw rate, roll attitude and roll rate. The complete state 


4-5 


required to describe the lateral modes of motion includes the observed component, 
sideslip angle. In the current DFCS. sideslip angle is obtained from guidance-and- 
navigation estimates or sensors with no additional filtering applied. As in the case 
of angle -of- attack, additional filtering will be applied to sideslip angle in subsequent 

versions of the DFCS- 

The lateral estimator also includes a "yaw-rate washout filter. This filter 
has a similar function to the pitch- rate compensation filter discussed aboi/e. The 
yaw-rate washout filter acts as a high- pass filter on yaw rate so that the steady- 
state component of yaw rate which is induced in a coordinated turn will not be fed 

back to the lateral control law. 
b) Im plementation 


Definitions 

i|i G * 4* 

(yaw attitude error) 


r-, - r 

(yaw rate error) 

Reduced lateral 

G 

state vector: 

* 0 - * 

(roll attitude error) 


a 

1 a° 

(roll rate error) 


( 4 , , r, 6 and p are the actual values of yaw attitude, yaw rate, roll 
attitude and roll rate and + G . r c . * G and p c are the desired values iron, guidance 
of yaw attitude, yaw rate, roll attitude and roll rate.) 

Reduced stat^ measurement: 



r in* ^in and p in are ^ n P uts to tlie * r *~ ter Update Routine,) 
Computations 

Sequential State Vector Correction (at Fast Sampling Rate) 
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is the state estimate following the i th correction where 
i -- 1, 2. 3, or 4. 


are the gain vectors which are multiplied by the measure 
ment residuals. 


Sideslip Angle Update: 

0=0 (no filtering applied to 3. 3 to be filtered at fast rate in 

p out in . rurr'c \ 

subsequent versions of Drtb.) 

Washout Filter Computation : 

, „ r + r (r is computed in the Filter Pushdown Routine) 

r wash = 1 • ' 5 out + r sum ' sum 

Computation of Roll Angle About Velocity Vector 

= 0 

v out V in 

* is required by the control laws. In the current version of DFCS, *v in 
is obtained'di rectly from the input interface routine and no transformation is required 
to obtain * v t . In subsequent versions of DFCS. <* Vout will be obtained from a 
transform^ applied to ^ and * out or the filter computations for the lateral state 

will be revised to supply * Vou t directl y as an output of the State estimator * 

C. Attitude Control Laws 

1. Longitudinal Control 
a) Engineering Description 

Longitudinal control is provided by tandem elevon deflection and by reaction- 
control jets. The tandem elevon deflection will be described as an "elevator” deflec- 
tion in the following discussion. 
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The objective of the longitudinal control system is to control vehicle pitch 
attitude to respond to commanded angle -of -attack from guidance. This is accom- 
plished through integral control of angle-of-attack. A commanded elevator deflection 
is computed in the DFCS control routine as a linear combination of three longitudina. 
state components: the integral-of-angle-of-attack error, angle-of-attack error and 
pitch- rate error. As discussed earlier, the control gains which multiply these 
errors have been pre- determined and scheduled as a function of Mach numbet and 
dynamic treasure. The control gains are updated by the DFCS during flight by 
table look-up and interpolation of stored data points (see discussion in Parameter 
Estimation and Update Section). Limiting is applied to each of the three inputs to 
the elevator control law to prevent extreme elevator excursions and excessive 
angular rates in response to large changes in the angle-of-attack command. 

The jet control law for pitch control is designed to aid the elevator control 
law in preventing large rate buildup. Control jets are programmed to fire only 
when the magnitude of pitch rate exceeds 5 /sec. 

b) Implementation - Aerodynamic Contro l Law (Elevator Control) 

Computation of integral 


A / o = 0. 5 


out 


n-1 


) T i 


"out / * F 
n 


out 


out 


n-1 


Limit ing of / or 

If / a > +4 deg-sec — » / or = +4 deg-sec 
If / o < -4 deg-sec -* / <* * -4 deg-sec 


Limiting of a Input to Control Law 
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lim 
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J * O , w + 


out* 


out 

' 3 -° 0 « a ’out< +3 -°° 
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i i miting of a Inout to Control Uw 


- t . 0° sec, q out = * 5 * 0 ' sec 


' . 0°/ sec < a. < 5.0 'sec 


+5.0° ‘sec. q out £ 5.0 sec 


Commanded Elevator Deflection 


6e C = * K l' /dr * K 2 * a ° ut lira * ^ ’ ^lur, + 6 * C 


(where 6e r 


is constant during transition and is initialized to the value 


of 6e_ existing at the conclusion of the entry phase) 

Speed Brake (Rudder Flare) Deflection 

6sb c =45° (Constant during transition phase) 


Notes: 


1) The indices n and n- 1 refer to the present and previous sample 
times respectively. 

2) Tp is the sample period (0.1 sec). 

3) In subsequent versions of the DFCS. determination of fa will 
be done by t e state estimator. 

4) Quantities not subscripted refer to present sampling lime. 
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In shaded band, no pitch jet firing is commanded. 

if . n > +o°/ 5 ec iet firings are commanded for a sample period 
4 out 

( 0.1 sec) to produce a negative pitch torque. 

If - a < - 5 0 sec, iet firings are commanded for a sample period 
4 out 

( 0.1 sec) to produce a positive pitch torque. 

2 . Lateral Control 
a) Engineering Description 

Lateral control is provided by differential elevon deflection, rudder deflection 
and reaction-control-jet firings. The differential elevon deflection will be described 
as an "aileron" deflection in the following discussion. 

The objective of the lateral control system is to control vehicle roll attitude 
about the velocity vector in response to commanded roll angle from guidance, and 
to maintain tight 3 ideslip control. Roll control is achieved by integral control of 
roll angle about the velocity vector. Yaw -sideslip control is achieved by feeding 
sides Up angle and yaw rate to the lateral aerodynamic controllers and to the jet 
control law. Both the commanded aileron deflection and commanded rudder deflec- 
tion are computed as a linear combination of five lateral control state components, 
yaw- rate error, integral -of- roll- angle error, roll -angle error, roll- rate error and 
sideslip-angle error. The lateral control gains are updated in the same way as the 
longitudinal control gains. Limiting is applied to the integral of roll angle and to the 
roll angle input to prevent extreme actuator excursions and excessive angular 
rates in response to large changes in the roll angle commands. 

The jet control law for lateral control consists of a channel to maintain tight 
sideslip angle control through the use of a sideslip-angle/ yaw -rate phase plane and 
a channel for roll control which employs a roll -angle/ roll -rate phase plane. 

A further aid to lateral-directional stability is obtained by a fixed rudder 
flare (speed brake) deflection of 45°. The 45° rudder flare deflection provides the 
maximum restoring yaw torque in response to sideslip error. 

b) Implementation - Aerodynamic Control Law (Aileron and Ru dder Control) 
Computation of ♦ Integral 

o.s(* + ® v ) T f 

out n-l ° ut n 
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- 6 


out 


out 


n-1 

/ 0 - / <*> . + A J ^ 

J n n- 1 

Limiting of J ^ 

If J0.' +3 deg-sec — */ d = + 3 deg-sec 
If / 0 < -3 deg-sec— *f d = -3 deg-sec 

Zeroing of J d for Small Roll Angle Errors 

If 1/41 >0.5 deg-sec and Id I < 0.125°- 

out 


/ d = o 


Limiting of d Input to Control Law 


v lim 


-3-0°. ^v out ® -3*0 C 


vout* ' 3 . 0 °< d vout <3.0° 
o 

-v < II 

v out 


+ 3.0°, >3.0 


C ommanded Aileron Deflection 


+ K • 3 + 6 

-r rs. - 0 ( 


«>C* K 9 ' r »» S h + K 10 • +K n ■ * v , inl * K 12' P out + 13’ out 


Commanded Rudder Deflection 


Sr C = - K l4' r wash' K 15' ; ^' K 16‘ 6 v 


- K,„ * P, 


lim 


1 7 ”out ! 8 out 


K ifi * 


+ 5r^ 


Notes: 

1) The indices n and n-1 refer to the present and previous sample times 
respectively. 

2) Tp is the sample period (0.1 sec). 

3) t n subsequent versions of the DFCS. determination of I * n will be done 
by the state estimator. 

4) Quantities not subscripted refer to present sampling time. 
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c) Implementation - Jet Control Law for Yaw and Roll Channels 


E R 



Notes : 

1) For the yaw channel, E = and ER = ‘ r wash * 

2) For the roll channel, E = -*v ou t and ER = 'Pvout 

{where py .is roll rate about the velocity vector and is presently 
being computed in the control section as: 

P v ^ * p out coso M + r wash sintt M }' 
out 

3) A given E, ER input will correspond to a point lying in one of the 12 
zones indicated in the above figure. 
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4) in zones 1. 3, 4. 5. « and 8 a rale change is commanded to drive the 
rate to the target rate indicated. For example. *n zone l, u 

ER ^2°/ sec fhen a rate change of 7. 73°/“ sec will be commanded. 

5) In the "coast'' zones 2. 7 and U a rate cnangc or zero is commanded. 

6) In zones r > and 11 a rate change is commanded to drive the rate to zero. 

7) In zones 10 and 12 a rate change is commanded to drive phase plane , 
location to the indicated sloped line bounding zone M. Tnese hnes save 
a slope of -1 sec-1 and intersect the axes at symmetrical dead baud. 

In the roll channel the deadbands are set at tl ( andtl ,o C . In 
yaw channel, they are set at td.2V>and *0.2 - see. 

3> The c >m.manded rate changes are meaificd before being passed to the 
jet select logic as follows: 

denote roll channeled rate change by 2 J r0 ]i anc * denote yaw 

_ ft- - t . J eiknrtiVA hv A . 


channeled rate change by 


vaw 


If lA ^roll !< °* 25 ° /sec 


,Aw roll = ° 


If 


yaw 


! > 0.03°/ sec 


If 


roll 


>0 

vaw 


» 0 

yaw 

^-roll 1 ° 


(The last constraint Is introduced to prevent a roll-jet firing in those 
cases in which it will produce a 3 ideslip-angle- rate opposite to that 
produced by the desired yaw jet firing.) 


D. Filter Pushdown 

a) Engineering Description 

The Filter Pushdown Routine of the DFCS for transition contains the extrapo 
lations of the current longitudinal and lateral state estimates to the next (fast) 
sample time. At that time, the extrapolated estimates are used in the sequential 
update equations in the Filter Update Routine described earlier. 

b) Implementation - Longitudinal Estimator 

^ out “ ^out + q out 

q out = q out + G 21 i G 22 q out + G 23 G 24 U out } 

where: 

Sq . c 

G s | 9B— Tp (where q dp is dynamic 

21 l yy * {pressure 

G 22 3 C M q 

G 23* C M a 
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(7H * mach nund>er) 


■5? 

% 

i 

T 

jf 


4“ 


I 

* 

i 



i 


t 

3 


"1 


i 




me 


M 


G 


23. 


2 4 


The coefficient si G M G 24 Are updated 
interpolation between -bored data points (see 
Update Section h« low ) . 


by the [UVS during flight bv linear 
di^ussion in Parameter F.stinnt'on and 


c ) Implementation ~ Lateral Ratimator 


* out " v out ^ 1 2 * out 1 F 

r ku 1 r out * C 2s| C, 26 'out * G 27 p out + °28 ^out I 
*out * *out * P out r F 

p mit ' ’’out * °29|°30 p oul * °31 r out * °32 8 out | 
■Vm * ->•” r out ’ O' 68 r »ash 


where: 


xx 


25 

. — - 
a 

I “ 

i 2 > 


* XX 

zz 

xz 


= ±-l 

f 

r* 

+ i*z, 

’26 

' 2V\ 

k 0 

i r 

•xx 


- h | 

ic 

+ L- 

*27 

~ 2V \ 

K 

! xx 


— Sq . hT « 

2 . ! d p p 


G«q “ C "t t ^ | 

28 n £ * xx 0 

1- 


xz 


zz 


StU„ b T 1 


' 29_ (1 1 -,T,-dp“-F 

u xx A zz xz' 


°30*^(%^ Cn p) 

° 31 *^( C « r ^ C . r ) 


Notea: 


32 


C, + ^C 

a ! zz n a 


1) The factors 1. 75 and 0. 98 are constant values in the current DFCS. 

2) The factor 0. 98 yields a time constant of 5 sec. 
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The coefficient* 0„ - 0 „ arc updated by the DECS during (ugh, by liner 
interpolation between etjred date points t,.r discussion in Parameter Estimation an 

Update Section below). 

^ Parameter Estimation and Update 
a) Engineering Description 

The computation of filter gains and the associated update of the state-covariance 
matrices used for filter gain computation are made in the Parameter Estimate and 
Update Routine of the DFCS. In addition, the coefficients required for state extra- 
polation and the control gains are determined through a table look-up procedure. 
Three value* of each control gain and extrapolation coefficient are stored in the 
table. The three values of each coefficient are the values associated with three 
different Mach numbers. The value of a coefficient for any arbitrary value of Mach 
number is determined through simple linear interpolation with respect to M.cn 

number. 

Each of the stored control gains contains an implicit multiplicative factor of 
dynamic pressure. When a value of control gain is extracted from the table through 
interpolation, it must then be divided by dynamic pressure to obtain the desired 

value. 

All of the parameter estimation and update equations are executed at the 
alow sampling period (2 seconds)* 

b) Implementation - Longitudinal Filter Ga i n and Covariance Estimation 
Definitions 


°36 °37 


G 37 G 38 


= Covariance of longitudinal reduced state, 0, q. 


a ? * Noise variance of 9 measurement. 

V 


s Noise variance of q measure ment . 


G 49 G 50 


G 50 G 51 


Process noise covariance matrix 
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Filter Gain and Covariance Update Computations 


1 

a 

n i 

1 

G ;tfi 

G„ 

* "h 

G '<7. 

L -J 




_ 


m - 


^ / / “ 

G tr, c, :i7 


1-G 1 ” 


°!r, C, :<7 

G 37 C, '58 

tm 


f 

i 

O 

ro 

L ^ 


G 37 C, 38 


- “ 



C 3 

1 

G 37 


‘ !" 


.°4. 

< G 38 + °q 1 

G 38 


G 36 

G 37 


1 

I 

a 

Vi 

1 


C 38 

G 37 

G 37 

G 38 


0 

1- G 4_ 


G 37 

G 38 


Covariance Extrapolation 


G 36 C 37 


H 

1 


’ C 36 g :<7 


1 o 


G 49 G 50 










G 37 G 38 


0 1 


G 37 G 38 


Jr 


G 50 G 51_ 


Note: 


In the implementation of the current DFCS. the above equates £r filter 
gain and covariance update are bypassed. The values o y . q ' 49 • 50 

and G 51 have all been set to zero. The gains Gj. G 2 , G 3 and G 4 are set 
to 1. 0, 0, 1 respectively. 


c) Implementation - Lateral Filter Gain and Covari ance Estimation 


Definitions 


| G 39 

O 

O 

G 41 

G 42 

o 

O 

G 43 

G 44 

G 45 

C 41 

G 44 

G 46 

G 47 

C 42 

G 45 

G 47 

G 4®« 


covariance of lateral reduced state, 

r, ^ 
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a ~ s Noise variance of b rneasui ement 



a ~ = Noise variance of r measurement 
r 

a “ 5 Noise variance of 0 measurement 

0 



Fitter Gam and Covariance Update Computations 
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Lateral Covariance Extrapolation Computations 
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c v 4 

T r . 

1 

0 

n 

1 

0 


F 




0 

0 

: 

t f 


0 

0 

1 

0 






0 

l' 

Tp 

1 

0 

0 

0 

1 




1* 





- 


_ 





Note : 

In the implementation of the current DFCS, the above equations for filter 
gain and covariance update are bypassed. The values of the measurement 
variances , * 2 , and ® p 2 and the components of PN have all been 
set to zero. The^ains G 5 -G 2() have been set to the following values: 

1. 0 , 0 . I, 0 , 1 , 0 . 0 . 0 , 0 . 1 . 0 . 0 . 0 . 0, 1 . 


d) Implementation - Confutation of Cont rol Gains and State Extrapolatio n 

Data is sto-.H for 3 longitudinal gains. 10 lateral gains and 10 coefficients 
used for longitudinal and lateral state extrapolation. Additional capacity for 5 
coefficient* is incorporated in the table in the current DFCS. The data is stored 
for three Mach numbers (0.9. 2.0. 5.0) and the parameters are determined by 
linear interpolation. This process is illustred by the figure below. 



Mach 

Number 


The numerical values currently stored in the DFCS are shown in Table 4-1. The 28 
coefficients and gains are denoted by V j — V 28’ 
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Table 4-1: Stored Data for Computation of Control Gains 
and Extrapolation Coefficients 


Coefficient 


0 * 9 


Mach Number 

2.0 


5 . 0 


V 1 

810.1 

V 2 

631.6 


425.7 

V 4 

0 

V 5 

0 

V * 

0 

V 7 

37.4 

v 8 

51.4 

V 9 

307.2 


150.0 

10 

v ll 

- 425. 9 

V 12 

81.2 

V 13 

- 8.0 

V 14 

- 51.7 

V 15 

- 18. 1 

V 16 

65.3 

V 17 

0 

V 18 

0 


0 

v 19 


V ^ 

0 

v 20 


V 21 

0 

V ^ 

0 

22 


V 23 

0 

V 24 

0 


0 

25 


V 26 

0 

V 27 

0 

V « 

0 

28 



* 180.0 

1300.0 

2287.4 

0 

0 

0 

37.4 

51.4 
307.2 
300.0 

0 

81.2 

- 8.0 

- 51.7 

- 18.1 

65.3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


89.0 

530.3 

327.1 
0 

0 

0 

30.58 

54.37 

325.1 

450.0 

400.0 
86.0 
- 8.44 

- 54.7 

- 19.2 

69.1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


1) V , V, , V R are not required for the current DFCS and are stored as 

4 9 5 6 

zeros. 

2 ) v.,- Vgg are data for state extrapolation equations. These equations 
are not implemented in the current DFCS and the data have been stored 


as zeros. 
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Eouations for Control Gains and Extra polation Coefficients Co m 
rJTm fWrait of linear Interpolation 


longitudinal gains : 


i 


K l = V l /q dp ; 

K, - V., q <|p ; 

? 

2) 

lateral gains: 


1 

* 


K -) ‘ V 7 q dp ' 

K 1P ’ V 8 q «lp * 

I 


K 12- V 10 /q dp 

' K n 'n q dp 

i 

! 

: 


K n’ V .2 , '<dp 

• K 1S 1 V l.V q dp 

5 

j 


K u* v i 5 'V 

* K 18 * V 16 q dp 


V V q c 


K 1 1 V 9'* q clp 


longitudinal extrapolation coefficients : 
°22 =V 17 ’ G 23 = V 18 * G 24 ' V 19 

lateral extrapolation coefficient?: 


G 21 ” V 26 T F 


G 26 = V 21 * 

G 27 ~ V 21 * 

G 28 = V 22 ’ 

G 25 = V 27 


G 30 ” V 23 * 

G 3i : ^ 24 * 

G 32 = ^2fS * 

G 29 " ^28 
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F. Block Diagrams of Traoaition Control Channel 

Figures 4-2, 4-3. anci 4-4 present a detailed control- engineer's description 
of the three control channels (Pitch, roll and yaw) by means of functional block 
diagrams. The computation of the control gains shown in these figures is described 
in the preceding subsection. 



Fig. 4-2 Transition Phase Pitch Loop 













Fig. 4-3 Transition Phase Roll Loop 














a> 

-o cn 

c c 



a 

c 

o 

-I 


fO rti 

cn J5 -o 

<D 

fO <D CT 

sl| 

a? o 
c c 

a> — -3 

-C w 
cn 

■o w <2 

<d o ~ 

C 

e a> 

05 is 
1 —iS) 

a aj r 
£ uo g* 

-a e 

' <= _C 


0 

c 

c 

15 

0L5 


E 

k- 

V. 

C“ 

*^__r 

3 

to 

lO 

O 

<d 


40 

CD 

TO 

aj 

-C 

O- 

3 


JC 

to 

c 

0 

W- 

CD 

Ic 

ro 

O 


05 

c 

10 

T3 

a> 



T3 

a? 

«— 4 

^05 

sz 

., _ 

CD 

»— 

O 

to 








SECTION 5 


ENTRY PHASE 

by 

J. Edwin Jones, Peter S. Weissman, Greg L. Zacharia» 


A. Introduction 

The DFCS design presented in this section is intended to provide stability -axis 
attitude control for the vehicle during entry, where this mission phase begins at 
approximately 400,000 ft altitude and ends at Mach 5. The control logic makes use 
of both the ACPS and the aerodynamic control surface torques, blended together in 
such a manner as to minimize ACPS fuel, while providing the capability for simul- 
taneous parallel use when necessary. 

The basic design features include gain scheduling and control sub-mode 
switching (both as functions oi selected trajectory parameters) to account for vari- 
ations in the vehicle dynamics and to allow for reasonably constant closed-loop 
transient response characteristics. In addition, except for very low dynamic 
pressure operation, maximum use is made of the aerodynamic control surfaces to 
minimize ACPS fuel. Finally, the controller synthesis, and thus the design pre- 
sented here, makes advantageous use of the vehicle's aerodynamic torques resulting 
from trim attitude deviations, allowing the vehicle's rigid body aerodynamics to aid 
in attitude maneuvers. 

Figure 5-1 provides an overview of the entry phase control-logic structure. 
Currently, measurements are supplied directly to the Control Routine by the Input 
Interface Routine. State estimation for the entry phase is in development. 

Longitudinal control of angle-of-attack and pitch rate is provided by elevator 
(tandem elevon) deflections and pitch ACPS jet firings. The elevator logic uses 
conventional rate and position feedback (with feedforward trim integration) through 
scheduled gains to maintain relatively constant closed-loop pitch dynamics, while 
the ACPS logic uses phase-plane control to maintain a desired limit-cycle operation. 

Lateral control of sideslip and roll angle is provided by aileron (differential 
elevon) deflections and roll and yaw ACPS jet firings. The aileron logic uses rate 
and position feedback of the lateral variables (with feedforward trim integration for 
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center- of- gravity offset compensation* to provide attitude and rate control in side- 
slip. The desired sideslip attitude, in turn, is determined to provide the roll ^roue 
necessary for roll control. Thr ACPS logic uses phase-plane control to maintain 
a desired limit- cycle operation when there is low dynamic pressure, and sideslip 
rate damping to maintain lateral stability when there is nigh dynamic pressure. 

In both lateral and longitudinal channels, blending -f oturol by the surtace.-. 
and the ACPS is provided through a logic which determines ACPS inhibition on the 
basis of the magnitude of the current surface comm irv!. 


B. Longitudinal Control 


1. Engineering Description 

Longitudinal control uses both elevator (tandem elevens) and pitch ACPS jets 
to maintain the commanded angle of -attack. There are two control sub-mod* s 
for the elevator, and one for the pitch jets. 

Early in the entry, when the dynamic pressure is low. the elevator is used 
strictly for trim, so as tc minimise pitch disturbanc* torques and thus Tuel usage. 
The trim control law is an open- loop quadratic function of the commanded angle - 
of-attack only. During this flight regime, the pitch ACPS jets, controlled by an 
angle of- attack/ pitch- rate phase- plane logic, provide a titude and rate control m 
response to guidance commands. There is no jet inhibition under this flight 
condition. 


When the dynamic pressure gets high enough, the elevator is used for both 
trim and transfer t control. Trim i. maintained through the use of an integrator 
acting upon the aagle- of- attack error; the initialization value of the integrator is 
calculated from rhe last trim value of the previous sub- mode (when dynamic pressure 
was low). Elevator control of transient errors is maintained through linear feed- 
back of angle -of- attack and pitch rate, with appropriate gains. During this flight 
regime, pitch jet control identical to that of the previous regime may be 
utilized together with the elevator, or control with the jets may be inhibited. This 
is done by the blending logic which is essentially a two-sided deadbnnded relay 


with hysteresis. The commanded elevator deflection is compared with preset 
fractions of its maximum and minimum values; should the command be "small”, 
it is as siimed that the elevator has sufficient control authority, and thus the pitch 
jets are ;iot required. Conversely, should the elevator command be large . the 
pitch jet control is not inhibited, but is allowed to assist the elevator in controlling 
the vehicle. This blended control continues throughout the rest of the entry, so that 


the pitch jets are always available for control assistance. 


The logic provides for a gradual turn-on of the elevator control so as to avoid 
a switching transient when the second control sub- mode is entered. 
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a) Elevator Control 


Figure 5-2a is a block diagram of the elevator control loops. When the dynamic 
pressure is less than q d . T . only the open-loop trim setting shown in the lower 
right-hand corner of the diagram is sent to the elevator. The function used to deter- 
mine the trim value is shown in Fig. a -2b. 


This parabolic curve is obtained as a lea^t- squares curie tit to the mono nt 
balance relation between angle-oi-attack and elevator deflection 1 the hypersonic 
regime. Note that the function is independent of both Mach number and dynamic 
pressure. 


When the dynamic pressure reaches q ( jpj. closed- loop control of the angle- 
of- attack is provided by conventional position (a land rate (q) feedback through a 
network which attempts to maintain a constant transient response throughout vhe 
flight envelope. Specifically, the desired closed- loop transfer function is chosen 
to be a well- damped second-order system, or 



desired 


s ♦ 2 


d U 'd 


S ♦ u!. 


(5-1) 


where <f d . * d > are specified for acceptable performance (k a is not explicitly 
specified;. The gains shown in the diagram may be separated into one of three 
types: <1> prespecified by the desired closed-loop characteristics; (2) dynamic 
pressure dependent <l/q dp >. or (3) scheduled gains which are functions of the aero- 
dynamic flight regime <e.g. . ( j). f Trim control if. accomplished in this sub- 
mode by a conventional clamped trim integrator acting on the angle-of- attack 
error to provide a steady- state trim elevator command in the absence of transient 
arigle-of-attack errors 

As can be seen from Fig. 5-2a and 5-2c, in order to allow a gradual blending 
in of active elevator control, the elevator command limits are functions of dynamic 
pressure which gradua’ly open up from the open-loop trim setting to the physical 
deflection limits. 

b) Pitch ACPS Control 

Shown in Fig. 5-3 is the ACPS phase-plane switch logic representation used 
in control of the pitch jets. For longitudinal control, the phase-plane coordinates 
are angle- of- attack error and pitch rate, or: 

(< 1 * f 2* * ^ (5-2a) 


See Subsection D for the numerical values of parameters. 
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The phase* plane is separated into three regions (by the switch curves), in which 
the commanded jet firing is either positive (u - +U>. negative (u - -U>. or «ro 
(u * 0). Thus for the particular application of pitch jet control 

(u. V) - <u y . l- y > (3-2b) 

whore u and V symbolize the jet command and the available acceleration level. 

y y 

respectively. 

The equations for the switch curves themselves are based or. the fuel-time 
optimal solution to the double integrator control problem and are modified by 
deaoband incorporation to ensure practical limit-cycle convergence (see typical 
trajectory of Fig. 5-3). The switch curves are defined by: 

>* : f i “ - S! ( 2 ) 

<5-3a) 

r ± : *1 “ ***2 + 2U f 2 * 

where, in this case, the deadbands are pitch attitude deadbands, or: 

(*!.«*> * 4 « 2 > (5 ' 3b> 


and ? is a fuel-time weighting constant greater than unity given for pit^h control 
by: 


<j * a 

or 


(5-3c) 


The interface (or blending) logic between commanded pitch jet firings and the 
commanded elevator deflections is shown in Fig. 5-4. The approach taken is to 
inhibit firing when there is "sufficient" pitch control acceleration from the elevator. 
This measure of sufficiency is obtained by comparing commanded elevator ( <Se c shown 
shown in Fig. 5-2) with threshold values to determine whether or not "excessive" 
elevator deflection is being called for. The hysteresis path included in the blending 
interface avoids chatter in elevator and ACPS activity due to cross-coupling effects. 
Note that this design includes an additional switch to preclude jet inhibition when 
the dynamic pressure is low and the elevator is used solely for trim. 

2. Implementation 

Coding implementation for longitudinal control during the entry phase is 
summarized in flow-chart form in Fig. 5-5. 


5-7 






Fig. 5-5 Logic Flow for the Longitudinal Control Channel in Entry Phase 
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c. 


Lateral Control 


1 . Engineering Description 

Lateral control uses both aileron (differential elevons) and roll and yaw ACPS 
jets to maintain small sideslip angles and follow commanded roll angle maneuvers. 
There are two control sub- modes for aileron, one for the >aw jets, and hree 


the roll jets. 


Early in the entry, when the dynamic pressure is low. the aileron command 
is 7 ero. and no lateral trim is attempted. During this regime, the yaw jets, 
controlled by a roll angle/yaw -rate phase-plane logic, provide attitude and rate 
of the roll ancle in response to guidance commands. The roll jets are 


controlled by one of two logics, depending on the yaw jet activity. When there are 
no yaw commands, the roll jets, controlled by a sideslip/ sideslip- rate phase- 
plane logic, maintain small sideslip attitudes and rates. When the yaw jets are 
commanded to maneuver the vehicle, the roll jets are controlled by a sideslip/ 
sideslip-rate phase-plane logic, which is similar to the first sub-mode except 
that compensation for yaw jet torques is included. Again, the purpose of this logic 
is to maintain small sideslip attitudes and rates. There is no jet inhibition under 

this flight condition. 


Later in the entry, when the dynamic pressure is between two specified 
values, the aileron command is still zero. The yaw jets are controlled by the same 
logic as in the previous regime. The roll jets are controlled by one of two logics, 
again depending on yaw jet activity. When there are no yaw jet commands, the 
roll j<*ts are commanded by a deadtanded relay logic which simply provides rate 
damping in sideslip. When the yaw jets are commanded to maneuver the vehicle, 
the roll jets are commanded by the same corresponding roll jet logic of the previous 
regime. There is no jet inhibition under this flight condition. 


Still later in the entry, when the dynamic pressure is high enough, the aileron 
is used for both trim and transient control. Trim is maintained through the use 
of an integrator acting upon the roll angle attitude error; the initialization value 
of the integrator is zero. Aileron control of transient errors is maintained through 
linear feedback of roll angle and sideslip attitudes and attitude rates with appropriate 
gains. During this regime, roll and yaw jet control identical to that of the previous 
regime may be utilized together with the aileron, or control with the jets may be 
inhibited. This is done with the same type of blending logic used longitudinally. 

Thus, with "small" aileron commands, lateral jet conti -a is inhibited. This blended 
control continues throughout the rest of the entry, so that the yaw and roll jets are 
always available for control assistance. As in the longitudinal case, there is a 
gradual turn-on of the allowed aerodynamic control authority. 
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a) Aileron Control 


Shown in Fig. 5-Ha is the block diagram tor aileron control^ 
control of vehicle roll angle is provided by conventional position (* 


C iost*d-loop 

a ) and rate 
w v 


(p, r) feedback of the lateral state variables through a 
maintain a constant response throughout the trajectory, 
closed-loop transfer function tor roll angle response ;s 


network which attempts to 
Specifically, the desired 
chosen to have four specified 


comp Lev poles, or: 




k ® (s ~ ~ tPl 


desired 


<s 2 + 2 


1 l 


S + a. ") (s 4 * - 2 “2 * 




(5-4) 


where (C-, ) are chosen for acceptable transient response (k 0 is not explicitly 

specified and is a function of the airframe dynamics.). This response is obtained, 
in the absence of rudder effectiveness, by taking advantage of sideslip mistrims to 
roil the vehicle. 

As with the longitudinal controller, the gains of Fig. 5-6 may be separated 
into one of three types: (a) prespecified by the desired closed-loop characteristics 
(e. g. , aj ); (b) dynamic pressure dependent; or (c) scheduled gains which are 
functions of the aerodynamic flight regime (e. g. , £ 3 ). For convenient reference, 
the gains of the first category are defined in terms of desired response character- 
istics with the following equation set: 


- , 2 , 2 

a l * u ' i “2 

a 2 = 2 (C p 2 + ^ 2 u, l >/(u) l U ' 2 ) 
a 3 = 2 '^ 1*1 + s 2 ^ 2 * 

a 4 = " l 2 + *2 + 4 4 1*2 


(5-5) 


where, as noted above, the (<?*., uj.) are chosen for acceptable roll response. 

As shown in the block diagram, lateral trim is accomplished through the use 
of a trim invegrator in the roll angle error channel. This allows the vehicle to trim 
to a non-zero sideslip angle (by commanding a non-zero trim aileron deflection) in 
the presence of body-axis roll disturbance torques due to lateral displacements in 
the vehicle*s center -of- gravity. 

The gradual expansion of the aileron command limits as dynamic pressure increase 
is shown in Fig- 5-6a and 5-6b. 


5-12 



















b) Yaw and Roll A CPS Control 

Yaw jet A CPS control synthesis is based on double integrator modeling of the 
roll angle dynamics, so that the switch logic presented for pitch ACPS control may 
be used. Roll jet logic consists of three models, one based on double integrator 
modeling of the sideslip dynamics, a second based on oscillator modeling, and a third 
which accounts for control axis cross-coupling, 

A double integrator model of the roll angle dynamics, combined with the 
assumption of turn coordination provided by the roll jets, allows for the use of 
the ACPS pitch jet control logic, with tne appropriate redefinition of control para- 
meters. Specifically, yaw jet control is defined by the phase plane switch logic of 
Fig. 5-3 and Eq (5-3a), where the coordinates are given by: 

(e., e 2 ) = t-® y sina r , ?) 
lim 

and tiie switch curve parameters are given by: 

(u, u> = <u 2 .iy 

* * 

The blending inter face logic is discussed telow. 

Roll jet control of vehicle sideslip is accomplished by the use of three control 
inodes. Two of these use a double integrator model for sideslip; the remaining one 
uses an oscillator model. 

The purpose of the first mode is to provide turn coordination when the yaw 
jets are firing (u * 0) . This is accomplished by the phase plane logic of Fig. 5-7, 
which is a generalized (3, 3) plane. The figure illustrates a typical state trajectory 
and the resultant limit cycle (with switch curve overshoots). The equations for the 
switch curve u are: 


(5-6a) 

(5-6b> 

<5-6c) 

<5-6d) 


= ^ [<*“>-* ] ^ S3>0 

r<& 


IT) 


SIH) 2 _ 1 


if S3 < 3 


m 


n : S3 


(5-7) 


here - is apptoMima*ril from the both rates: 

-* - p sin - r ro- j.j, 

hr switch curve parameters of Eq (3-7) are defined by 
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Fig- 5-7 Roll Jet Control when u z * C 



S - sg'\ (u .) 


\n\ 


( 2 V., * 


(5-9) 
(Cent > 


2 1 


where is* defined below in Eq{>-L0). 


The purpose of the «. < otul mode is to provide direct control over sideslip 
when the yaw jets are inactive in y - 0) an.l when the dynamic pressure iv low. 
With double integrator modeling, the phase-plane logic of Fig. 5-3 and Eq a- >a) 
may be used, where the coordinates arc given by. 


(tj.tj) » (P.a> 


<J-iOa) 


and the switch curve parameters are given by: 


(u. to * <’i x . U x sin» T ) 

(5-10b) 


(5-10c) 

7 = 

(5-10d) 


It is convenient to label this control mode (used at low q dp and when u„ , - 0 ) as 
attitude hold, since both rate and attitude error are driven toward zero. 

The third control mode takes advantage of the vehicle's dynamic stability 
in sideslip, which becomes significant once the dynamic wcsBure becomes suffi- 
ciently high. V/h-'t is required is simply rate damping of the sideslip angular 
velocity as shown in Fig. 5-8. As above, it should be noted that k is approximated 
by the body rates of Eq (5-8). 



Fig. 5-8 Roll Jet Control when u z = 0 
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The * utitrol modes described above, defining lateral control with yaw and 
roll ,ot. . are summarized in bloc k diagram form in Fig. &-«• Several points 
Should J. noted. First, the mode select block which defines the mil jet control 

mode to be used is defined by: 


” z ' q dp v q dp;; " (5 - ll) 

U Z 3 °* q dp iq dpn ‘ n 

where the modes are coordination, attitude hold, and rate damping, respectively. 

The three roll control blocks correspond with the above inodes. A second point to 
note is the gating by the blending parameter X a . *hich l* defined by logic that 
is identical to that of the pitch channel. Finally, it may be noted that the logic dis- 
allows ACPS inhibition during aileron tnm-only (null aileron command) operation. 

2. Implementation 

Coding implementation for lateral control during the entry phase is presented 
in flow-chart form in Fig. 5-10. 

D. Parameter Estimation 

All parameters used in control computations for the entry phase are either 
fixed or scheduled on the basis of Mach number, 771 . or mean (trim) angle-of-attack 
a . Table 5-1 tabulates all entry- phase fixed control parameters. The scheduled 
parameters, control gains ^ through £ g , arc each plotted as a function of the 
appropriate scheduling parameter in Fig. 5*11. 

It an initialization pass through the entry control logic, several parameters 
which are functions of pad- loaded values arc calculated and stored for future use. 
The Load Manipulation Routine for the entry phase is shown in Fig. 5-12. 
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Fig. 5-9. Yaw and Roll ACPS Control Block Diagram 
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Fie. 5-10 Logic Flow for + he Lateral Control Channel w Entry Phase 
(Cont) 
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10 Logic Flow for the Lateral Control Channel in Entry Phase 
(Cont) 
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Table 5-1 Fixed Parameters U3ed in Entry Phase Control Computations 


VARIABLE NAME SYMBOL 


AX 

u < 

2.8 

deg/ sec 

AY 

u 

y 

0.628 

deg/ sec 

AZ 

U z 

1.15 

deg / sec 

DAI 

*1 

0.75 

deg 

DA2 

*qf2 

1.0 

deg 

DB1 

*8 1 

0.2 

deg 

DB2 

*82 

0.4 

deg 

DR1 

Si 

0.75 

deg 

DR2 

S 2 

1.0 

deg 

DYNP1 

q dpl 

2 

lb / ft 2 

*> 

DYNP2 

q dp2 

20 

lb /ft 

DYNP3 

q dp3 

5 

lb /ft 

DP1STR 

q dpl* 

10 

Ib/fT 

9 

DP2STR 

q d P 2* 

30 

lb /ft“ 

-1 

KALFT 

K o T 

0 . 2 

sec 

-1 

KPHIT 

K <t>T 

0.01 

sec 

R1 

r l 

5 

deg 

R2 

r 2 

-45 

deg 

R3 

r 3 

+ 15 

deg 

R8 

r 8 

25 

deg 

R9 

r 9 

20 

deg 
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able 5-1 Fixed Parameters Used in Entry 


Phase Control Computations (Cont) 


VARIABLE NAME SYM POL 


VALUE 


SIGMA A 

a 

Of 

SIGMA B 

°3 

SIGMA R 

% 

TF 


W1 

“l 

W2 

u 

2 

WD 


ZETA1 

’i 

ZETA2 


ZETAD 



D 


2 


5 

— 

0.1 

sec 

0.7 

-1 

sec 

0.4 

-1 

sec 

0.7 

-1 

sec 

1.4 

— 

0.7 

— 

0 7 

— 
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Fig. 5-12 Entry Load Manipulation 






SECTION 6 


ON-ORBIT PHASE 


by 


Donald W. Keene, Edward T. Kubiak (MSC/EG2) 


The DFCS (Digital Flight Control System) will control the attitude and the 
translation of the Space Shuttle during the various on-orbit activities: coasting 

flight, orbital maneuvers using the Orbital Maneuvering System engines or the 
Attitude Control Propulsion System (ACPS) engines, docking, etc. This capa- 
bility has not yet been incorporated into the DFCS program. However, a pre- 
liminary design for switching logic to control attitude error using the ACPS has 
been developed. To best perform the variety of on-orbit control functions 
demanded of the DFCS, two interrelated switching logics are proposed rather 
than a single compromise switching logic. The two algorithms are defined with 
respect to function and design features in the remainder of this section. 

Nominal Logic 

The function of the nominal logic is to provide attitude and rate control 
whenever disturbances can be neglected. This includes all operations except the 
docking maneuver (for which there is separate logic). The logic is *hown in Fig. 6 
and has the following features: 

1) Special coast zone helps to assure minimum impulse limit cycles; 

2) Placement of switching curves should minimize attitude overshoot 
due to damping of maneuver rates; and 

3) Rate limit region assures acceptahle rate levels for all operations. 

The nominal switching logic is partitioned into seven switching regions as 
shown in Fig. 6-2. In regions as I. II. III. IV and VI the control acceleration is 
negative. In regions V and VII the control acceleration is positive. 
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A summary of the desired control actions for each region is shown in Table 6*1. 
"Hysteresis" logic is included in the switching code to assure completion of the 
commands issued in regions HI. IV. VI and VII. This eliminates the possibility 
of switching line chatter. Typical phase point trajectories illustrating the hysteretic 
nature of the coast zone are shown in Fig. 6-3. A table of symbols used in the 
switching logic is given in Appendix A. Flowcharts for the nominal logic are given 
in Appendix B. 

B. Disturbance Logic 

The second switching logic is a special logic for attitude control in a 
disturbance environment. Typically, the disturbance would be due to gravity 
gradient torques and/or thrust misalignment torques during either OMS (Orbital 
Maneuvering System) or ACPS translational burns. Two separate logics are used 
depending on the size of the position and rate errors. For large errors, where the 
position error exceeds PL or the rate error exceeds RL. the logic shown in Fig. 6-4 
is used. This logic drives the state into the small error control region oy 
maintaining a constant rate. The small error switching logic is shown in Fig. 6 -5. 
The important advantages of this logic are: 

1) One can fix the limit cycle based on any one of the following considerations 


td T 2 
8 


by choosing the appropriate deadband. 

a) To fix limit cycle period, T: 



fc> To fix limit cycle maximum rate, 
* 


X 

max 


v Kl (*m» J 

s '* DB * 2 I u dI(I L cI TdI) 


c) To fix limit cycle width, simply set X DB to the desired value. 

2) By properly biasing the right and left deadband values, DR and DL, the 
time average attitude error can be set to zero. For example, let 
be the desired attitude variation. Then for a positive constant disturbing 
acceleration, the deadband values are: 

DR = X DB -X 

DL = 
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Fig. 6-4 Large Error Disturbance Logic 

If the position error magnitude exceeds or the rate error 
magnitude exceeds then the state is forced into the small 
error control region by rate limiting control. 
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Fig. 6-5 Disturbance Logic Switching Regions 
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Where X is the '’centroid” of the limit cycle and is derived in Appendix D 



3) this logic cures two of the problems observed m gravity gradient limit 
cycles during Apollo lunar orbit attitude control. 

a) Poor convergence due to straight-line switching 

b) Convergence to only minimum impulse firings when longer 
firings would be more efficient. 

The desired control actions for each region of the disturbance logic phase 
plane are summarized in Table 6-2. As in the nominal switching logic, hysteresis 
is used to assure completion of commands issued in regions I and 11 in order to 
avoid chattering about the switch lines . Typical state point trajectories are 
illustrated in Fig. 6-6. The symbols are defined in Appendix A. The flowcharts 
for the disturbance logic are given in Appendix C. 

C. General 

The following comments apply to each switching logic: 

1) The sample period. T. will be constant and is tentatively 0. 1 seconds. 

2) To eliminate the possibility of chattering along switching lines (e.g. , due to 
inexact control acceleration modeling) '’hysteresis' 1 regions are defined 

in the coast region. When th~ st. te is found to be in a hysteresis region, 
an additional check is made to c ermine whether control action was 
required in the preceding iteration (i. e. RFLG>0). If it was, controx 
is continued until the state is i.u longer in the hysteresis region. In this 
manner thruster firing times are updated each cycle and discrete error 
spikes {e.g., those due to CDU transient on apply) cause spurious firings 
only in the cycle in which they occur. 

3) The switching logic parameters are selectable on a per axis basis (e.g. , 
deadbands). 
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i 

i 

! 4) As presently forseen the disturbance logic will be used: 

! a ) when there are significant slowly changing disturbances (i. e. , 

| gravity gradient or aerodynamic) during long term attitude hold 

f b) during significant non-gimbaled engine thrusting maneuvers 

( caused by thrust misalignment) 

! \s these disturbances may not significantly effect all axes, the disturbance 

- logic ts selectable on a per axis basis. 

* 

i 
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APPENDIX 6A 


Symbol 


X„ 


u. 


D 


*N’D 
LRL, URL 

RDB 

RL. PL 

DB, DB’, DB" 



DR, DL. DL’ 



I 


ARL 

X D 

MAX 

Q. P 

R4F 


TABLE OF SYMBOLS 


Definition 

Units 

Position control variable 

degrees 

Rate control variable 

o , 

/ sec 

Control acceleration 
Disturbance acceleration 

o . 2 

/sec 

o / 2 

/sec 

Deadband factors 

degrees 

Lower and upper rate limits 

O / 

/ sec 

Rate deadband 


Rate and position limit 

°/sec, degree: 

Deadband values used 
nominal logic 

degrees 

Scale factors 

non- dim 

Disturbance logic position 
deadband 

degrees 

Deadband values used in 
disturbance logic 

Scaled accelerations 

degrees 

o. 2 

/sec 

Centroid of disturbance 
logic limit cycle 

degrees 

Index to denote x, y or z 
control axis 

non-dim 

Upper and lower rate limit 
average 

o . 

/ sec 

Desired rate change 

o t 

I sec 

Rate change corresponding to 
full cycle jet-on-time 

o. 

i sec 

Intermediate computational terms 

degrees 

Flag which allows portion of 
Region IV between URL and 
RDB to be a coast zone 

non-dim 
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APPENDIX 6A 


TABLE OF SYMBOLS (Cont) 


Symbol 


AR 


RFLG 


RAV 


K 


1 


K 2 


RR 


H MAX 


Definition 


Units 


Width of rate limiting coast ! sec 

zone in large error disturbance 

logic 

Flag denoting regions of the non-dim 

phase plane from which a 
command is issued (one per axis) 

Same as ARL 1 sec 

Flag used to denote the sign of non-dim 

the attitude rate.(Since the 
switching logic is symmetric 
about the origin the computa- 
tion of the desired rate change, 

X D* is performed by trans- 
forming the phase point so 
that it lies above the Xj axis. 

The correct sign is affixed to 
Xp) by multiplying by before 
exiting. ) 

Flag used to denote the sign of non -dim 

the attitude rate on the previous 

cycle (used to prevent chattering 

about the zero rate line when 

the phase point lies in the hyster- 

etic coast zone). 

_ o , 

Rate error used when phase / sec 

plane command is zero. (Used by 
jet selection logic is selecting 
appropriate roll, pitch jets. ) 

Rate limit used is coast zone to / sec 

govern jet selection. 
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FLOWCHART OF NOMINAL SWITCHING LOGIC 
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APPENDIX SC 

FLOWCHART OF DISTURBANCE SWITCHING LOGIC 




R&j 


DISTURBANCE LOGIC 
INITIALIZATION: 


INPUT: |U c |. RL, PL. C D< A D> C L> X DB . A I 
ALL QUANTITIES OF DIMENSION THREE 


ENTER FOR EACH AXIS 


INPUTS from ot^er 
ROUTINES 

X 1 

• POSITION ERROR 

x 2 

- RATE ERROR 

U, 

* DISTURBANCE 

d 

ACCELERATION 


COMPUTATIONS 

©K R . iu d l/|u c | ©DL' ■ DL- 

©5 • -3-(‘^) ©*» ’ 3|v d 


^)dr * X DB - x 
©DL - C L X 


2|u d I 
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APPENDIX 6D 


AVERAGE ATTITUDE ERROR DURING CONSTANT 
DISTURBING TORQUE LIMIT CYCLE 


A limit cycle of the type shown below is assumed. 


X 



DB 


The state is driven by a disturbance acceleration, Ujy along the trajectory 
from A to B. At point B a control acceleration, U c , is appliedreturning the state 
to point A. The problem is to find the average attitude error. X. in terms of the 
attitude travel. X„ D . Knowing X, a switching logic can be defined which translates 
the limit cycle along the Xj-axes such that the resulting X is zero. 

X will be found by determining the average attitude error for each of the 
trajectories OB and BO' (called Xj. and X 2 > and then "weighting” them according 
to time to calculate the net average, X. 

Determination of X^: 

Xj is defined as follows: 

*1 

x i -t $ X* 2 * ®' J) 

1 o 
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tj can be calculated from the Xj ordinate of B 2 • Note that: 


Bx , ’ 5 U D ‘! 2 


(D-2) 


B is found by first rearranging Xj * fj(X^) and X^ - fgfXg/ such that. 


X 2 " gjtXj). X^ * g^Xj) 
Then, at point B, gj(X^) * ggfXj), and, 

ft' - U, 


and from Eq (D-3) 


B x, ■ x db( C t- J - ) 


. . -/ 2X DB <"c - V 
1 t c tD 


Using t t in Eq (D-l) 


t . i * / 0 «-M 

4 T *DB y J 


(D-3) 


(D-4) 


(D-5) 


<D-6) 


Determination of Xjt 

The Hme average attitude error for the trajectory segment BO' can be found 
by setting t * 0, at point B. Then, 


* 2 =-r \ 

2 O 


or. 


X 2 * X o + “TT + i ' U c H 2 


(D-7) 


X is defined in Eq (D-4) and X can be determired from, 
o ° 


X o* U 


. . - + / 2X DB U D <U c- V 

j d‘i >/ tj- 


(D-8) 


To find t„ one can use the following relation, 
6 


X o + (U D- t 'c ), 2"° 
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• 1.1 


«MF«' 






From which it follows. 


t 


2 


_/ 2X db u d 


(D-9) 


Substituting the results of Ea (D-4), (D-8) and (D-9) into Eq (D-7), yields: 

,(> - ) <d -‘°' 


X 2 = X DB! 


Determination of X 


Finally, the total average attitude error. XT can be found from: 


= * 1*1 + * 2*2 
X “TT — 


(D-ll) 


‘1 *2 

Substituting the results of Eq (D-5), (D-6), (D-9), (D-10) into Eq (D-ll), yields: 


TL 



(0-12) 
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